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Abstract 



The transition 7*(<Zi)7*(<72) — * 71-0 (p) is studied within the QCD sum rule framework. As 
Oh: a first step, we analyze the kinematic situation when both photon virtualities are spacelike 

Oh. and large. We construct a QCD sum rule for i ? 7 * 7 * vr o (qf, g|) an d show that, in the asymptotic 

limit |<7i|, — ► °o, it reproduces the leading-order pQCD result. Then we study the limit 
\q\ | — > 0, in which one of the photons is (almost) real. We develop a factorization procedure 
for the infrared singularities ln(g 2 ), 1/qf , l/<?i , etc., emerging in this limit. The infrared- 
sensitive contributions are absorbed in this approach by bilocal correlators, which can be 
also interpreted as the distribution amplitudes for (almost) real photon. Under explicitly 
formulated assumptions concerning the form of these amplitudes, we obtain a QCD sum 
rule for F 7 * 7 .„-o(g^ = 0, g| = ~Q 2 ) an d study its Q 2 -dependence. In contrast to pQCD, we 
make no assumptions about the shape of the pion distribution amplitude (p n (x). Our results 
agree with the Brodsky-Lepage proposal that the (5 2 -dependence of this form factor is given 
by an interpolation between its Q 2 = value fixed by the axial anomaly and 1/Q 2 pQCD 
behaviour for large Q 2 , provided that one interpolates to a value close to that dictated by 
the asymptotic form <£>" s (x) = 6f n x(l — x) of the pion distribution amplitude. We interpret 
this as an evidence that <p n (x) is rather close to the asymptotic form. 
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1 Introduction 



The studies of the transition form factor for the process 7*7* — > vr° when two virtual photons 
7* produce a neutral pion is apparently the cleanest case for testing QCD predictions for elastic 
processes. In contrast to the pion electromagnetic form factor, the perturbative QCD subprocess 
7* (91) + 7*(<?2) —* q(xp) + q(xp) appears at the zeroth order in the QCD coupling constant a s , 
and the asymptotically leading term has no suppression. The relevant diagram resembles the 
handbag diagram for the forward virtual Compton amplitude used in the studies of deep inelastic 
scattering. This gives good reasons to expect that perturbative QCD for this process may work at 
accessible values of spacelike photon virtualities q\ = —q 2 ,q 2 = — Q 2 . In the lowest order, pQCD 
predicts that [|TJ 



where ^p- K (x) is the pion distribution amplitude and x, x = 1 — x are the fractions of the pion 
light-cone momentum carried by the quarks. In the region where both photon virtualities are 
large: q 2 ~ Q 2 >lGeV 2 , the pQCD predicts the overall 1/Q 2 fall-off of the form factor, which 
differs from the naive vector meson dominance expectation F 7 * 7 * 7r o(g 2 , Q 2 ) ~ l/q 2 Q 2 ~ 1/Q 4 - 
Thus, establishing the 1/Q 2 power law in this region is a crucial test of pQCD for this process. 
The study of F 7 * 7 * 7r o(g 2 , Q 2 ) over a wide range of the ratio q 2 /Q 2 of two large photon virtualities 
can then provide a nontrivial information about the shape of tp n (x). 

However, experimentally most favourable situation is when one of the photons is real q 2 = or 
almost real. For experiments on e + e~-machines having one small virtuality strongly increases the 
cross section. It was also proposed to use real photons to study the n° production in 'je collisions 
0. On fixed-target accelerators, like CEBAF, one can attempt to study the F J * J * 7V o(q 2 = 0,Q 2 ) 
form factor through 7*7r* — > 7 processes || with a photon in the final state. In the q 2 — > 
limit (with Q 2 large), the nonperturbative information about the pion is accumulated in the same 
integral / of <f 7r (x)/x that appears in the asymptotic pQCD expression for the one-gluon exchange 
contribution for the pion electromagnetic form factor. Hence, information extracted from the 
studies of F 1 * 7 * n o(q 2 = 0, Q 2 ) can be used to settle the bounds on the pQCD hard contribution to 
the pion EM form factor. 

Since the zeroth moment of the pion distribution function is normalized by the matrix element 
of the axial current (i.e., by the pion decay constant f n ), the value of I is sensitive to the shape 
of the pion distribution amplitude <p n (x). Two most popular choices are the asymptotic form 



(pl s {x) = 6f v x(l-x) [§-§ and the CZ model ^ z (x) = 30/^(1 - x)(l - 2x) 2 0. Using ^ z (x) 



increases the integral by an extra factor of 5/3 compared to the value based on (p°^(x). This 
observation can be used for an experimental discrimination between these two models. In fact, 
both the CELLO data || and preliminary high-Q 2 CLEO data || seem to favour the leading- 
order pQCD prediction with the normalization corresponding to a rather narrow distribution 
amplitude close to (p® s (x). To perform a detailed comparison of the (future) data with theoretical 
predictions, one should take into account the pQCD radiative corrections. These include the 




(1.1) 
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one-loop contributions to the hard scattering amplitude JTT], [12], which, decreasing the leading- 
order result by about 20%, still leave a sizable gap between the predictions based on two models 
mentioned above. One should also add the terms generated by two-loop evolution of the pion 
distribution amplitude [ jT5] , |H| , |Tj|. Originally, these corrections were found to be tiny A 
recent progress |16| in understanding the structure of the two-loop evolution suggests that the size 
of these corrections is somewhat larger. However, the numerical analysis of the two-loop evolution 
presented in ref.[17] does not indicate appreciable changes for the integral over the distribution 
amplitude. Hence, there are good chances that the controversial subject of the shape of <Pn{%) 
may soon be settled experimentally. 

Within the pQCD approach, the pion distribution amplitude ip. K {x) is a phenomenological 
model function whose shape should be taken either from experiment (this was not possible so far) 
or calculated in some nonperturbative approach, e.g., using QCD sum rules. However, applications 
of the QCD sum rules to nonlocal hadronic characteristics (functions), like distribution amplitudes 
(p(x) are much more involved than those for the simpler classic cases [TE| of hadronic masses and 
decay widths. The trickiest problem is that the underlying operator product expansion (OPE) for 
the relevant correlators has a slow convergence because some terms are parametrically enhanced 
by powers of N for the x N moment of the distribution amplitude 0. For this reason, one needs a 
very detailed information about the nonperturbative QCD vacuum to get a reliable sum rule. Such 
information is not available yet, and results for <p n (x) have a strong model dependence [p^j. In 



the present paper, instead of starting with the pQCD approach and taking (p n (x) from QCD sum 
rules, we calculate F^* y * n o(q 2 = 0, Q 2 ) directly from a QCD sum rule for the three-point function. 
A serious problem for such an attempt is that one of the photons has a small virtuality and the 
relevant three-point amplitude is sensitive to nonperturbative long-distance QCD dynamics. For 
this reason, as an intermediate step, we construct a QCD sum rule for a simpler kinematical 
situation when both photon virtualities are large. To take the q 2 —>■ limit, we perform additional 



factorization using the methods developed in refs.[20|-|27 



The paper is organized as follows. In Section 2 we introduce our basic object: the correlator 
of two vector and one axial current, discuss its behaviour in different kinematical situations and 
interrelation between QCD sum rules and pQCD approaches. The situation when both photon 
virtualities are large is considered in Section 3. We construct there the QCD sum rule valid 
in this "large-g 2 " kinematics and analyze its structure and some limiting cases. In Section 4, 
we outline specific problems one faces trying to use the operator product expansion in the limit 
when one of the photon virtualities is small (small-g 2 kinematics). The basic features of mass 
singularities, which appear in this limit, are illustrated using a scalar model as a toy example. 
The methods developed there are then used in Section 5 to construct a factorization procedure 
of long- and short-distance contributions in the realistic QCD situation. We show that the long- 
distance contributions in this case have the structure of bilocal correlators^ which are considered 
in Section 6. In particular, we discuss there the continuum and p-meson contributions into the 
bilocal correlators. In Section 7, we study the contact terms [20] which appear in some bilocals. 



The bilocal correlators were originally introduced in a similar context by Balitsky ES 
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In Section 8, we collect together the contributions calculated in preceding sections and write down 
the QCD sum rule for the -F 7 » 7 * 7r o(g 2 , Q 2 ) form factor valid in small-g 2 kinematics. We analyze 
the q 2 — ► limit and present our results for F 7 . 7 .„-o(g 2 = 0,Q 2 ). In the concluding section we 
summarize our findings. Some technical details of our calculations can be found in the Appendix. 



2 Form factor 7*7* — > tt° and three-point function 
2.1 Definitions 

The form factor i 7 7 * 7 .^ 7r o (g 2 , q 2 ) of the 7*7* — > tt° transition is determined by the matrix element 
4tt J (tt,P \T{J^X) J v (0)} \0}e'^ x d 4 X = V2ie m2 F rr ^o (q^q 2 ) (2.1) 
where J M is the electromagnetic current of the light quarks (divided by the electron charge): 

J » = (J«7m m - g<#V*) , (2-2) 
\tt,P) is a one-pion state with the 4-momentum p and we use throughout the convention e^apq^ = 

^-^vaq 2 i ^iMvaflQl Q.2 = ^lLvq\q 2 i &tc. 

To study the form factor, we should construct first a formalism in which the pion would emerge 
as a QCD bound state in the qq system. A possible way is to start with a three-point correlation 
function 

^071,92) =27ri J^\T{jl{Y)J^X)JM}\^ iqiX e^ Y d A Xd A Y , (2.3) 

(cf. pU| ) where p = q± + qi- In addition to the two EM currents present in eq. ( |2.1|) , the correlator 
( |2.3| ) contains also the axial current 

f a = (ui5laU - d>y 5 >y a d) . (2.4) 

The latter has the necessary property that its projection onto the neutral pion state is non-zero. 
In fact, this projection is proportional to the tt" — ► \iv decay constant /„. « 130.7 MeV: 

(0\f a (0)\7T°,p) = -iV2Up a . (2.5) 

The three-point correlator ( |2.3| ) has a richer Lorentz structure than the original amplitude 
(|2.1|), and not all the tensor structures it contains are relevant to our study of F 7 * 7 *_ >7r c(g 2 , q 2 ). 
Incorporating the Lorentz invariance properties of the three-point function and Bose symmetry 
for the virtual photons, one can write the amplitude T ailiV as 

F al Mu(qi, 92) = P a e OTl «^ (V, ql, qt) + r^^^Ai (p 2 , g 2 , qfj 

+ [tanq iq2 qiv — £avq 1 q 2 q2ij] ^2 (p 2 ', <2l, qt) + [ e a Wl q 2 q2u ~ ^auq^Ql^] -^3 (p 2 , Ql, <?2 
+ {^a miq2 qiv + e Q i/gig2<?2At] ^2 (p 2 , 9l> qf) + [ e a f iq 1 q 2 q2v + ^auqxq^lfj] -A 3 (p 2 > 9l> #2, 



ijy2 q^.2 

y -T± + y^ 5 + (pr)A 6 



+ e 



a [lisp 
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where p = q% + q 2 , r = q% — q 2 . The invariant amplitudes T, A have the following symmetry 
properties: Ti (p 2 , q\, qf) = Ti (p 2 , q\, q\) , A (p 2 , ql qf) = ~A (p 2 , ql <??)• 

Utilizing the fact that, in the four- dimensional space-time, there is no antisymmetric tensor of 
rank 5, i.e., 

^-afiwy 9Se ^fivyS 9ae ^wySa Qfie ^^/Safi Que ^SafjLU 9"f£ > (^-6) 

and using the conditions q\T ailv = q^T^u = imposed by the EM current conservation, we get 
finally the expansion in terms of four invariant amplitudes: 



T oifj.ui.qi, q 2 ) 



[ Pa F + r a A] 



+ [q 

+ e 



-ctfiur 



auqiq2i 

„2 



p 2 + r 2 p (pr) 



F + [q 2u e 



auqiq2 

(pr) 



A 



F + l±^A 



According to eqs. (|2.1| ), Q2.5|) , the three-point amplitude T a ^ u (qi, q 2 ) has a pole for p 2 



(2.7) 



2 

ml: 



T ^u(Qi, 92) — ~Y^Z_ 2 'Po£livqxq< 2 F~ f *~ f *-+K° (?lJ I2) + ■ • ■ (2-^ 

p 



Thus, the tensor structure of the pion contribution is p Q e^„ gi(?2 . 111 eq.(|2.7p, it corresponds to the 
invariant amplitude F (q 2 , q 2 ,p 2 ). In other words, F (qf, q 2 ,p 2 ) has the pole l/(p 2 — rn 2 ), and the 
form factor F 7 * 7 *_ i . 7r o (qf, q 2 ) can be extracted from the residue of that pole. 

However, the problem is that the bound state poles are present only in full amplitudes, formally 
corresponding to the total sums over all orders of perturbation theory (in QCD one should not 
forget to add also nonperturbative contributions). Terms corresponding to any finite order do not 
have such poles. Fortunately, it is not always necessary to perform an explicit all-order summation 
(impossible in QCD) to extract information about a particular bound state. 



2.2 Implications of the axial anomaly 

Using the axial anomaly [0 relation for massless quarks 



e 2 



we obtain the constraint 

P 2 F - (qf + q 2 2 )F + (q\ - q 2 2 )(A + A) = -. (2.10) 

7T 

For real photons (q\ = q\ = 0) , this reduces to p 2 F — 1/ir, provided that F and (A+A) do not have 
1/q 2 or l/q 2 singularities, which is true since there are no massless qq bound states in the vector 
channel. Hence, F ~ 1/p 2 for real photons, i.e., F really has a pole corresponding to a massless 
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pion pl| , p2| . Furthermore, the anomaly relation ( |2.10| ) fixes the value of the F 1 *- y *^ n o(qf 7 q 2 ) 
form factor when both virtualities of the photons are zero : 



7 7 



^(0,0) 



(2.11; 



It should be emphasized that the anomaly relation does not require that F 7 » 7 »^ 7r o(0, Q 2 ) has this 
value for all Q 2 : eq.( [2.i0| ) can be satisfied even if F 7 . 7 .^„-o(g 2 , Q 2 ) has a nontrivial dependence on 
photon virtualities. 

Thus, the axial anomaly allows one to calculate exactly one particular combination of invari- 
ant amplitudes (|2.10|) for arbitrary values of the virtualities q\-,q\-,p 2 - In QCD, this is a rather 
exceptional situation. Normally, a reliable QCD calculation is possible only in the region of large 
virtualities where one can incorporate the asymptotic freedom property of the theory. 





Figure 1: a) Three-point correlation function, b) Structure of factorization in the limit 
Q 2 ,q 2 > \p 2 \. 



2.3 Factorizable contributions and perturbative QCD 

If both of the photon virtualities are large, the leading term of the l/g 2 -expansion of any diagram 
contributing to the amplitude F can be written in the factorized form (cf. |33], |34|): 



[2.12) 



where C(£, 77, q\, g 2 ) is the short-distance coefficient function and IT(£, 77,79) is the long-distance 
factor given by a particular term of the PT expansion for the correlator of the axial current with 
a composite operator 0(£, 77) ~ g(£) . . . q(i]): 



m,V,p) ~ / (0|T(O(e,7 / )j(r))|0)e i ^ 4 y. 



;2.i3) 
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Figure 2: a) Lowest-order pQCD term, b) Hard gluon exchange diagram. 
The correlator T],p) is formally given by a sum over all orders of PT, and it has a pole for 

2 2 



n(£,7?;p) 



p2 _ m 2 



<0|C(e,77)k°, P> . 



;2.i4) 



Comparing eqs. ( |2.14| ) and ( |2.8|) , we conclude that the transition form factor is given by 

F rr -^{<$,<$) = I C(^ V ,q 1 ,q 2 )(0\O^ V )\7r°,p) d^d 4 V . (2.15) 



Thus, the possibility to factorize short- and long-distance contributions allows one to get an explicit 
expression for the form factor of a bound state, the pion, without explicitly obtaining the pion 
pole from a summation or a bound-state formalism. 
Now, introducing the distribution amplitude <p n (x) 



/o 



where x = 1 — x, we obtain the hard scattering formula 



F rr ^o(ql, q\) = T(qx, q 2 ; xp, xp) <f n (x) dx , 
with T(qi, q 2 ; xp, xp) being the amplitude for the subprocess 7*7* — > qq. 

2.4 Perturbative QCD predictions 

Calculating the subprocess amplitude in the lowest order one gets the result |jj 

4tt f 1 (p v (x) 



[2.16) 



(2.17) 



F L °tf,Q 2 ) = — [ ^ [X) dx 
3 Jo xQ 2 + xq 2 



X18) 



which was already discussed in the Introduction. For defmiteness, from now on we will adhere 
to the convention that q 2 is the smaller of the two virtualities: q 2 < Q 2 . The leading-order 
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pQCD formula has a smooth limit when q 2 — > 0, predicting that the asymptotic behaviour of the 
77* — > 7T° form factor is [p: 

F^AQ 2 ) = f [ ^rfx + 0(1/Q 4 ). (2.19) 

The x n -moments of the pion distribution amplitude ^p n (x) are given by the matrix elements of the 
twist-2 composite operators g7s{7 1 'D Vl . . .D Un }q. Furthermore, the contributions to FjS,J?(Q 2 ) 
from the higher twist two-body operators ^^{D" 1 . . . D Un }q and q^^a^ D Vl . . .D Un q, which are 
potentially dangerous due to large magnitude 0(m 2 /(m u + m d )) of their matrix elements || [35| 
are proportional to the light quark masses, and their net input is suppressed by 0(m 2 /Q 2 ) factor 
for large Q 2 . Hence, one should not expect unusually large (1/Q 2 ) N corrections to the leading 
twist term (g7Tg ). 

In the first order in a s (Q 2 ), one can imagine a subprocess in which the real photon first disso- 
ciates into a quark-antiquark pair (the relevant amplitude can be called the photon distribution 
amplitude <p y (y)). The gg-pair then interacts with the virtual photon to produce the final state 
quark-antiquark pair converting eventually into the pion. This contribution is analogous to the 
pQCD hard scattering term for the pion form factor. However, in our case, such a subprocess is 
formally realized only at the next-to-leading twist level: 

F^° (Q 2 ) ~ a s (Q 2 ) f ^f^ dx + 0(1/Q 6 ). (2.20) 

A more careful inspection shows, however, that if (p^(y) ~ y for small y, the integral over y 
diverges. This means that this contribution is not factorizable in the usual sense and should be 
considered within the QCD sum rule context (see discussion below). We will see that this term is 
rather small, mainly because it is suppressed both by 0{a s /^) and 0{1/Q 2 ) factors compared to 
the leading-order term ( 2.19Q . 



One-loop perturbative QCD corrections to eq. (|2.19|) are known [10 - ||12|| and they are under 
control. For these reasons, the process 77* — > 7r° provides the cleanest test of pQCD for exclusive 
processes. 



2.5 Wave-function dependence and Brodsky-Lepage interpolation for- 
mula 

To fix the absolute magnitude of the leading-order pQCD prediction ( |2.19| ), one should know the 
value of the integral 

I= ±r?M dx , (2.21) 

which depends on the shape of the pion distribution amplitude <p n (x). In particular, using the 
asymptotic form|4]]-[|6| 

<p™(x) = 6Ux(l-x), (2.22) 
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one obtains I as = 3 resulting in the absolute prediction for the asymptotic behaviour of the 
F 11 *^o{Q 2 ) form factor [p]] 

F™ (Q 2 ) = ^ L + O(l/Q 4 ). (2.23) 



On the other hand, if one uses the CZ-amplitude, the value of the integral I increases: I cz = 5, 
and experimental data can, in principle, discriminate between these two possibilities. Of course, 
the asymptotic l/Q 2 -dependence cannot be a true behaviour of F 77 . 7r o(Q 2 ) for all Q 2 -values: it 
should be somehow modified in the low-Q 2 region to comply with the bound on F 11 * lx o(Q 2 ) in 
the Q 2 = limit imposed by the anomaly relation fl2.1ip . Brodsky and Lepage |36| proposed the 
interpolation formula 

F^\Q 2 ) = - r . 1 n ^ (2-24) 



which reproduces both the Q 2 = value (|2.11| ) and the high-Q 2 asymptotics Q2.23|) dictated by 



the asymptotic form of the distribution amplitude (|2.22| ). According to refs.|8], |(|, this formula 



agrees with experimental data. It also agrees with the results obtained in several quark model 
calculations |]3~?||-||40||. On the other hand, the curve based on the formula 



interpolating between the Q 2 = value and the CZ-normalized high-Q 2 prediction, is far from 
existing data points and quark model results. 

2.6 Pion wave function and QCD sum rules 

From the theoretical side, there are also doubts fl9| , f|l[] that QCD sum rules really require that 
the pion distribution function has the CZ shape. In particular, the QCD sum rule calculation of 
(Pn(x) at the middle point x = 1/2 performed in ref. [[4T| produced the value (p n (l/2) ~ 1-2/, 



7T- 



to be compared with (p°^{l/2) = 1.5 f n and <^ z (l/2) = 0. In ref. (TP]], it was pointed out that 
keeping in the OPE the lowest condensates only does not provide information necessary for a 
reliable determination of the pion distribution amplitude. This is especially clear if one writes the 
CZ sum rule directly for (pn(x) 

= ^a-e- S0/M2 )x(l-x) + ^l[5(x)+S(l-x)} 

+ ^^^- {U[6{x) + 6(1 ~ X)] + 2[6 ' {X) + 6 ' {1 ~ X)]} - (2 - 26) 



As emphasized in ||19H , it is the 5-function terms here which are crucial in generating a humpy 
form for (p^ z (x). Adding higher condensates, e.g., (qD 2 q), one would get even higher derivatives 
of S(x) and 5(1 — x). All the subseries of such singular terms can be treated as an expansion of 
some finite-width functions related to nonlocal condensates. The sum rule based on a model for 



9 



the nonlocal condensates consistent with the earlier estimates for (qD 2 q) |42| reduces the values of 
the lowest nontrivial moments of (p 7r (x) bringing them very close to those for the asymptotic form 
(|2.22| ). As a result, the model distribution amplitude constructed in |19[ gives / ~ 3. One can also 
try to construct a QCD sum rule directly for the integral /. However, it is clear that such a sum 
rule cannot be derived by a simple substitution of the original CZ sum rule into the integral ( |2.21| ), 
because of singularities generated by the S(x)/x and 5'(x)/x terms. The singularities disappear if 
one uses the nonlocal condensates in the sum rule for (^(x), and the resulting value for / is close 
to 3 f43j. A more radical way is to consider the sum rule for the original 77* — > it amplitude as 
a whole, without approximating it by the first term of the pQCD expansion. As we will see, the 
singularities mentioned above will appear then as the power (l/q 2 ) n infrared singularities in the 
relevant operator product expansion (OPE). However, these singularities are caused by a formal 
extension of the OPE formulas from the large q 2 region where they are valid, into the small-g 2 
region, where the OPE breaks down. In fact, the OPE should be modified in the small-g 2 region. 
Such a modification is equivalent to the regularization of the infrared singularities. Our goal in 
the present paper is to construct a QCD sum rule for the transition form factor valid in the region 
of small q 2 . 



2.7 Unfactorizable contributions and QCD sum rules 

If the photon virtualities are not very large, then it is normally impossible to factorize the p 2 - 
dependence of the diagrams contributing to the three-point amplitude F (p 2 , qf, q%) into separate 
factors one of which will produce the pion pole. We know, of course, that the full amplitude 
F (p 2 , q 2 , g|) must have the pion pole, but it is unclear how much a particular finite-order diagram 
contributes to such a pole. 

To display the pole structure of F (p 2 , qf, gf), it is convenient to use the dispersion relation for 
the three-point amplitude: 

F (p\ ql qf) = - r P(g,9l2, f 2) ds + "subtractions" . (2.27) 
v ' n Jo s — p 2 

The pion contribution to the spectral density is proportional to the F 7 * 7 *^ 7r o form factor: 

P ( s ) (Zi) Q2) = ^fir^i 3 — m^)F^» 7 *^o (qf, q 2 J + "higher states" , (2.28) 

but any approximation for F (p 2 , qf, q 2 ) obtained from a PT expansion, gives information not only 
about the pion, but also about the higher states. Hence, there are two problems. First, we should 
arrange a reliable PT expansion for F (p 2 , q 2 , g|). The only way is to take p 2 spacelike and large 
enough to produce a reasonably converging OPE. The second problem is that when p 2 is large, the 
pion contribution does not overwhelmingly dominate the dispersion integral and, to disentangle 
it, one should subtract the contribution due to higher states. Note, that though the contribution 
of a state located at s = m 2 is suppressed by the — p 2 / {m 2 — p 2 )-factor compared to that of a state 
at s = 0, the relative suppression of higher states disappears as — p 2 tends to infinity. 
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The higher states include A\ and higher pseudovector resonances which presumably become 
broader and broader, with their sum rapidly approaching the pQCD spectral density. So, the 
simplest model is to approximate all the higher states, including the A\, by the spectral density 
p PT (s, q\, gf) calculated in perturbation theory: 

P m ° d (s, ql qt) = irfJ(s)F rr ^o (ql q 2 ) + 9(s - s )p PT (s, ql g 2 ), (2.29) 

where the parameter s a is the effective threshold for the higher states. 

The smaller p 2 , the bigger relative contribution of the lowest state. Hence, the strategy is to 
take the smallest possible p 2 within the region where the 1/p 2 expansion is still legitimate. In 
fact, it is more convenient to use the faster decreasing exponential weight exp[— s/M 2 } instead of 



l/(s — p). This is achieved by applying to F(qf, q\,p > ) the SVZ-Borel transformation |T8| : 



B(p 2 - M 2 )F(ql q\,p 2 ) = $(g 2 , ql M 2 ) = jT e'^ P (s, ql q 2 2 ) ds. (2.30) 

Another merit of the SVZ-Borel transformation is that using it one gets a factorially improved 
OPE power series: l/{-p 2 ) N -> (1/M 2 ) N /(N - 1)!. 

In a sense, the QCD sum rules can be treated as a method of extracting information about 
the lowest state from the behaviour of F(qf,qlp 2 ) in the large-p 2 region. To construct a QCD 
sum rule, one should calculate the SVZ transform $(g 2 , ql M 2 ) as a power expansion in 1/M 2 for 
large M 2 . To this end, one should calculate first the three-point function Tip 2 , q 2 , Q 2 ) as a power 
expansion in 1/p 2 for large p 2 . However, a particular form of the 1/p 2 expansion depends on the 
interrelation between the values of the photon virtualities q 2 and Q 2 . 

The simplest case is when both virtualities are sufficiently large and similar in magnitude: 
Q 2 r*j q 2 no —p 2 > where p is a typical hadronic scale p 2 ~ 1 GeV 2 . This case will be referred 
to as the "large-g 2 " kinematics. Then all power-behaved (1/M 2 ) n contributions correspond to 
the situation when all the currents J M (X), J„(0) and j%(Y) are close to each other, i.e., all the 
intervals X 2 , Y 2 , (X — Y) 2 are small. 

A more complicated case is when q 2 is small q 2 <C p 2 , while Q 2 is still large: Q 2 > p 2 . 
In this case, to be referred to as "small-g 2 " kinematics, one should also take into account the 
configuration when the electromagnetic current J^{X) related to the g-photon is far away from 
two other currents, i.e., when there is a possibility of long-distance propagation in the g-channel. 
In the limit q 2 — > 0, such a propagation is responsible for the mass singularities [1/q 2 , hi(g 2 ), etc.) 
in the Feynman diagrams contributing to the three-point amplitude. 



3 QCD sum rules for the i^^o^ 2 , Q 2 ) form factor in large- 
q 2 kinematics 

Let us consider first the simpler case, when both Q 2 and q 2 are large. In this situation, it is 
sometimes convenient to introduce another set of variables: average virtuality Q 2 = (Q 2 + q 2 )/2 
and the asymmetry parameter u = (Q 2 — q 2 )/{Q 2 + q 2 ). 
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3.1 Lowest-order perturbative term 

The starting point of the operator expansion is the perturbative triangle graph (Fig.^a). Its 
contribution contains almost all types of the tensor structures present in eq. ( |2.7|) . Extracting the 
F part and using the Feynman parameterization we get: 

F PT {q\Q\p 2 ) = - f 1 — ^ - 5(1 - ^Xi) dx x dx 2 dx 3 . (3.1) 

7i Jo [q z xix 3 + Q z x 2 x 3 - p z X\X 2 \ ~y 

In this representation, it is straightforward to apply the SVZ-Borel transformation to get 

j>"V,Q 2 .M 2 ) = ^ t^{- qW TS™ 3 \ W-t*)**'!*'**- (3-2) 

Comparing this representation with the definition of the SVZ transform (|2.30|) , we can immediately 
write down the formula for the perturbative spectral density p PT (s,q 2 ,Q 2 ): 

p PT (s } q 2 , Q 2 ) = 2 f 1 5 ( s - q2xiX3 + Q2x2X A 5(l-j2^i)dx 1 dx 2 dx 3 . (3.3) 

Scaling the integration variables: x\ + x 2 = y, x 2 = xy, X\ — (1 — x)y = xy and taking trivial 
integrals over x 3 and y, we get 

p PT (s, q 2 , Q 2 ) = 2 f XH f 2 n + 2 l q Zu dx , (3.4) 
Jo [sxx + xQ z + xq z \- i 

or, in terms of Q 2 and uj: 

PTl 2 o f 1 XXQ 4 (1 + UJ(X - X)) 2 , . 

p^ 1 (s,q z ,Q z ) = 2/ _ 1 1 — dx. (3.5) 

Jo [sxx + Q 2 (l + oj(x — x))} 3 

It should be noted that the variable x in the integrals above is the light-cone fraction of the total 
pion momentum p carried by one of the quarks. 

A very simple result for p PT (s, q 2 , Q 2 ) appears when q 2 = 0: 



p PT (s,q 2 = 0,Q 2 ) = I -f W2 . (3.6) 



As Q 2 tends to zero, the spectral density p PT {s, q 2 = 0, Q 2 ) becomes narrower and higher convert- 
ing into S(s) in the Q 2 — > limit [ [3l|j : 



p PT (s,q 2 = 0,Q 2 = 0) = 5(s). (3.7) 

Thus, the perturbative triangle diagram "tells us" that, for massless quarks and Q 2 = 0, two 
photons can produce only a single massless pseudoscalar state, and there are no other states in 
the spectrum of the final hadrons. As Q 2 increases, the spectral function broadens, i.e., higher 
states can also be produced. 
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3.2 Gluon condensate corrections 

When the Borel parameter M 2 (or the probing virtuality — p 2 ) decreases, both perturbative (loga- 
rithmic or 0(a s ) ) and nonperturbative (power or 0(l/p 2 )) corrections come into play. As argued 
by SVZ fllSfl , the power corrections proportional to quark and gluon condensates: (0\ipi/j\Q) 2 , 
(01(7* (7" |0), etc., are much more important than the higher order perturbative corrections. In 
many cases, the latter can be safely neglected. 




d) e) f) g) 



Figure 3: a) Lowest-order perturbative contribution, b) — g) Gluon condensate corrections. 

The lowest-order diagrams proportional to the gluon condensate are shown in Fig.|3|&-g. They 
take into account the fact that propagating through the QCD vacuum, quarks interact with 
the nonperturbative gluon fluctuations which can be treated as a background field. The most 
straightforward way to calculate these diagrams is to take quark propagators in the coordinate 
representation using the Fock-Schwinger gauge z^A^z) = for the background field A p (z): 

+ 4^ y ' A '" G »* (0) " 19^SI (A ' 2y2 " ( Xi ') 2 )Gft(0)G, 3x (0). (3.8) 

Here A = X — Y and G a/ 3 = ^e al 3 pa G pa . Using this expression for each quark propagator of the 
original triangle diagram and retaining the 0(GG) terms, one obtains all the diagrams of Fig.|3]&-g. 
In particular, one can immediately see that contributions ^fl and |3|e vanish due to our choice of 
the coordinate origin at the axial current vertex (such a choice, in which the two photons are 
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a) b) c) 




d) e) f) 




g) h) i) 

Figure 4: Quark condensate corrections involving "soft" gluons. 

treated symmetrically, is more convenient in this calculation than that implied by our original 
definition ( p. 3D ). The remaining diagrams ^),c,f,g are easily calculated in the coordinate repre- 
sentation. After performing the Fourier transformation to the momentum space, we extract the 
tensor structure corresponding to the F form factor and then apply the SVZ-Borel transformation 
to the relevant contribution F^ GG ^(q 2 , Q 2 ,p 2 ). The final result for the sum of the 0(GG) diagrams 
reads: 

= —(—GG) ( -J- - J-) * (3.9) 
3.3 Quark condensate corrections 

For massless u- and <i-quarks, the quark condensate contribution starts with terms proportional 
to (0\qTqqTq\0) . Using the usual vacuum dominance hypothesis these can be reduced to 
(0|gg|0) 2 . There are two types of diagrams producing the (0|gg|0) 2 contributions. 
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First, there are diagrams shown in Fig.^a-z (sometimes called the soft gluon diagrams) cor- 
responding to local operators of qDDDq, qGDq and qDGq type, with the covariant derivatives 
DDD, and the GD, DG factors producing an extra qq term by equations of motion. 

The diagrams shown in Fig.|5]a-r (and corresponding to a hard gluon exchange) give the 
(0|gg|0) 2 structure directly. However, only the diagrams ^i-d contribute to the F form factor. 
The total 0((qq) 2 ) contribution to the SVZ-borelized amplitude is 



64n 2 a s (qqY 
243M 2 



M 4 



243M 2 



9 q 2 ' 
2Q 2 + Q\ 


9 

+ QV 


+ w) 


/ll -3u 2 


18 \ 

■ + 


l 


V Q 2 M A 


(1-cu 2 ) 2 



(3.10) 



3.4 QCD sum rule 

Combining now eqs.flggg), ( p3U|) , (gg), Q) and ( p30|) , we obtain the QCD sum rule for the F 
form factor valid in the region where both virtualities of the photons are large: 



nUF rrn o(q 2 ,Q 2 ) = 2 Tdse-^ 12 f 

Jo Jo 



xx (xQ^ 



xq 2 ) 2 



64 
243 



M 4 



9 7T 

Q 2 9 9 



[sxx + xQ 2 + xq 2 } 3 

1 1 1 

+ 



dx 



2M 2 Q 2 2M 2 q 2 Q 2 q 2 



+ 4 

g 4 2g 2 2Q 2 Q 4 



9 9 



QV QV 



(3-11) 



or, in terms of the variables Q 2 and u: 
F- 



+Hl/^LGG) 

9 7T 



- - — ►7r°(? 

1 



g 2 ) = 1 ( 2 f S °dse-^ M2 1' dx 
nf n [ Jo Jo 



xxQ 4 (l + oj(x — x)) 2 
sxx + Q 2 (l + u(x — x))] 3 
1 1 



1 



1 



,q 2 m 2 g 4 / 1 



or 



64 
243 



TT 3 az s (qq)' 



11 — 3c^ 2 18 

g 2 M 4 + Q£ 



2\2 



(3.12) 



In particular, taking the exactly symmetric kinematics, when uj = and q 2 = Q 2 , we obtain the 
sum rule 



-^7*7*— >7T° (Q )Q 



y 7r 



1 



1 



1 



dse~ s ' M2 
64 



1 xx Q 4 
dx 







[sxx + Q 



213 



(3.13) 



73 + 0T0 71 " a s <99>' 
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Q 2 M 2 Q 4 / 243 



18' 



g 2 Af 4 q 



derived in ref. 29 
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3.5 Q 2 — > oo limit and transition to perturbative QCD 

In the limit Q 2 — > oo with u; fixed, eq.( |3.12| ) reproduces the sum rule 



F rr ^o({l-uo)Q\{l + uj 



UQ 2 I t 1 U (l 



xx dx 



+ u>(x — x)) 



+ 



7T 



.a,--. 64vr 2 (11 - 3uj 2 ) . , 2 1 , o „ „, 



9M 2 (1-W 2 ) 7T 

obtained by Gorsky p4 |, who calculated the leading 1/Q 2 contribution only. 

To make a connection with the perturbative QCD approach, it is instructive to rewrite eq. ( 3.14j ) 

as 

4vr r 1 dx { 3M 2 



7 7 



^((i-^)g 2 ,(i + ^)Q 2 ) 



3f w Q 2 Jo (l + u(x-x)) \ 2tt< 



-(1-e 



-so/A/ 2 



I XX 



+^(^GG)[8(x) + 5(x)\ + ^vra s (gg) 2 (ll[5(x) + 5(x)] + 2[6'(x) + *'(*)])}. (3.15) 

Note, that the expression in large curly brackets coincides with the QCD sum rule ( |2.26| ) 
for f n <p n (x). Hence, in the large-Q 2 limit, the QCD sum rule ( |3.15| ) exactly reproduces the 
perturbative QCD result 

47T r 1 ip n (x) dx 



F rr ^o ({l-uj)Q 2 ,{l + uj)Q' 



3Q 2 Jo (l+co(x-x))' 



(3.16) 



4 Operator product expansion in small-g 2 kinematics 



In the nonsymmetric kinematics g 2 <C Q 2 ~ 1 GeV 2 , it is more convenient to use the "old r 
variables q 2 = —q\ and Q 2 = —q\ instead of uj and Q 2 . 



4.1 General features of the small-g 2 kinematics 

Incorporating a general analysis of the asymptotic behaviour of Feynman diagrams (see Appendix 
A and |3"3"|, one can show that, in the small-g 2 kinematics, there are two different regions of 




a) 




+ 




Figure 6: Configurations responsible for the power-behaved contributions. 
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Figure 7: Separation of SD- and S-contributions for perturbative diagrams and gluon condensate 
corrections. 



integration capable of producing a contribution to the correlator (see Fig.|6|), which behaves like an 
inverse power of p 2 . The first one is the purely short-distance region corresponding to a situation 
when all three currents are separated by short intervals, i.e., all the intervals X 2 , Y 2 , (X — Y) 2 
are small. The second region corresponds to another short distance regime, in which the vertex 
X related to the small virtuality photon is separated by a long interval from two other currents 
(this means that Y 2 is small, but X 2 and (X — Y) 2 are large). This is illustrated by Fig|| where 
Fig||z shows the full correlator J r a ^ tu ( ( li, whereas Figs|]6, c represent the two possibilities of 
getting the power-law contributions. 

In Fig.||c, the long distance contribution (the dashed blob) is given by a two-point (bilocal, cf. 
28| ) correlator 11(g) of the electromagnetic current J M and some composite operator O of quark 



and gluon fields, the latter represented by <g>. At low q 2 , the correlator 11(g) cannot be calculated 
in perturbation theory. The standard strategy is to model this nonperturbative object by the 
"first resonance plus continuum" ansatz, with the parameters of the spectrum determined from a 
QCD sum rule for such a correlator. In other words, one should calculate 11(g) at large q 2 using 
the OPE, then extract, in a standard way, the parameters of the model spectrum in the g-channel 
and, finally, use the model spectral density in a dispersion relation for 11(g) to obtain 11(g) at 
small q 2 . 

On the diagrammatic level, the factorization procedure should be applied for each contributing 
diagram: the lowest-order triangle graph (Fig.|7]a, b, c), gluon-condensate graphs (Fig. [Tjd, e, /), 
etc. In this way, each diagram is represented as the sum of its purely short- distance (SD) and 
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bilocal (B) parts, the latter obtained through the relevant contribution into the operator product 
expansion for J V (Y) j„(0). In fact, it is more convenient to define the SD-part of a diagram as the 
difference between the original unfactorized expression and its -B-part determined via the OPE. 
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The total SD-contribution is given by a sum of the SD-parts of all relevant diagrams: 




Figure 8: SD-contribution. 

Now, substituting the SD-terms by the "original minus 5-terms" , we obtain the expansion shown 
on Fig.|^. The dots there stand for the rest of the condensate diagrams and higher-order correc- 
tions. The first row in Fig|| corresponds to the usual OPE for the correlator ( |2.3| ) in the large-g 2 






+ 




p- 



Figure 9: Structure of factorization in the small-g 2 kinematics. 



kinematics 



q 2 ~ 



Q 2 



-p 



jj 2 ). The second row corresponds to additional terms which should 



-p 2 ~ jj 2 ). Note, that at 



be taken into account in the case of small-g 2 kinematics (q 2 Q 2 
large values of q 2 , these additional terms should die out to convert the modified OPE specific to 
the small-g 2 kinematics into the standard OPE for the large-g 2 kinematics (cf. |24], 56]). Another 
statement is that, at low q 2 , the "subtraction" terms (those in the parentheses in the second row 
of Fig||) have exactly the same singularities as the corresponding terms in the first row of Fig.||] 
and, therefore, the total expression is regular in this kinematic limit. 
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4.2 Logarithmic singularities 



According to eqs. ( |3. 11| ), the condensate terms contain l/q 2 , 1/g 4 , etc. singularities in the q 2 — > 
limit. On the other hand, the perturbative expression, though finite as q 2 — > 0, contains contribu- 
tions which are non-analytic at this point. To study the structure of the non-analytic terms, it is 
convenient to use the following method. First, let us introduce another set of integration variables 
in the expression for the Borel transform $ PT (g 2 , Q 2 , M 2 ) (]3T 



xi + x 3 = X, x 3 = yX, xx = (1 - y)X = yX, x 2 = 1 - A. 



This gives 



<^(g 2 ,Q 2 ,M 2 



xdXdye-^ x / M2 h~y Q2 'y M2 



(4.1) 



ttM 2 jo 

Performing a formal g 2 -expansion of the exponential in the integrand, we will get divergent in- 
tegrals for all the coefficients of the (q 2 ) n expansion, except for the lowest term. To get a more 
sensible result, we use a continuous version of the series expansion for the exponential, i.e., the 
Mellin representation: 



A J T(-J) dJ . 



(4.2) 



1 r-5+ica 

Z = TT-r / 

2m Js-ioo 

Now, the integral over A can be taken easily, and the next step is to calculate the J-integral 
by taking residues at J = 0, 1, . . .. This gives: 



<^ T (g 2 ,Q 2 ,M 2 



1 



7T 



M 2 Jo 



dye 



2 V/M 2 y I 1 , QV_ q*y/M* , 

+ M 2 + 



2 4 2 

M 2 M 4 



\M 2 



~ /^\>(n)( n + i; 



n 



1 ! 



(4.3) 



where ip(z) = T'(z)/T(z) and r(^) is the Euler Gamma function. 

The non-analytic terms q 2 ln q 2 and q 4 ln q 2 are typical examples of mass singularities ||46|| . They 
appear when a long-distance propagation of particles is possible. In our case, the mass singularities 
are related to the possibility to create, for massless quarks, a gg-pair by a zero-virtuality photon. 
To apply perturbative QCD methods, one should first factorize the contributions due to short and 
long distance dynamics, as outlined above (see Figs.|7|.|8|) . 



4.3 Scalar model 

To clarify the origin of the singularities and outline the method of their factorization, let us 
consider first the analogue of our form factor in a simple scalar theory g<f>hy. 

H<lu<l2) = J d 4 Xd*Y e-^ x e^ Y (0\T{j(X)j(Y)j(0)} |0>, (4.4) 

where j(X) =: <j)(X)<j)(X) : . The perturbative contribution and some of the power corrections 
are shown in Fig.[lO|. 
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b) 



c) 



d) 



Figure 10: Lowest terms of the OPE in a scalar toy model. 

Calculating the perturbative contribution (FigJTOa) in the a-representation or using the Feyn- 
man parameters, we get 



1 q 2 x\x 3 +Q 2 x 2 x R 

$ (a) = B(-p 2 ^ M 2 )F [a) = I dx 1 dx^9{x l +x z <l)- e 



dXdy 



o 



e -yq 2 \/M 2 \ e -yQ 2 /M 2 y 



X 2 X 3 



lo XyM 2 

Using the Mellin representation, as described above, we obtain the expression 



(4.5) 



(a) 



2tt 2 Jo 



dye 



— I 

yM 2 \ 



n=0 



yq 2 \ n il)(n + 1) 
M 2 



nl 



(4.6) 



in which the term containing the logarithmic singularity is explicitly displayed. 
The diagrams |i~0|&, c, d corresponding to power corrections give 



$ (6) = B(-p 2 M 2 



g 2 Q 2 



0, $ 



2 ) 



Q 2 M 2 ' 



g 2 M 2 



(4.7) 



Since the diagrams |10"|a, d contain terms singular in q 2 , it is necessary to perform an additional 
factorization of short- and long-distance contributions for these diagrams. 

To get the bilocal contribution for the amplitude ( |4.4j ), let us extract terms related to the 
simplest coefficient function proportional to the propagator S(Y) = i/A.ii 2 (Y 2 — iQ) (see Fig.||c): 



r 1 00 1 

T B = / d A Y e ipY -— Y — Y^ ... Y" n 
J tt 2 Y 2 ^ Q n\ 

x / d 4 Xe~^ x (O|T{ J (X)0(O)(^ 1 ...^J0(O)}|O). 
Here, the long-distance contribution is described by the two-point correlator 
DU^fa) = / d'Xe-^ x (O|T{jW0(O)Rx---flk}^(O)}|O> 

= (-«) n 9l M l • • • Qlftn ^(ft 2 ) + • • ■ , 

where the dots denote the terms containing g^^. 



(4i 



(4.9) 
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The large-p 2 behaviour of T B is governed by the leading light-cone singularity of the integrand. 
To display the F 2 -structure of the integrand in eq.(|4.9|), we reexpand the composite operator over 



the traceless operators possessing a definite twist p4| : 



[n/2] | / 9/ I "O / va2 \ ' 

= E ";r , .1 ~r ^ • • • ^ n - 2i H<K<9 2 yRi . . . d, n _ 2l ^ . ( 4 .io) 



„ l\{n-l + l)\ 



In ( |4.10|) , we used the standard notation {<9 Ml . . . <9 Mn } for a traceless group of indices /ii . . . fi n : 
grWWO{... w ... w ...} = 0. The leading 1/p 2 contribution in ( |4.8| ) in the p 2 —>■ oo limit comes from 
the lowest twist operators (t = 2 in this case). All other operators are accompanied by the 
factors ~ Y 2 , (V 2 ) 2 , etc., which cancel singularity of the propagator ~ 1/Y 2 and, hence, give no 
contribution in the large-p 2 limit. 

To get the B-regime contribution for the perturbative triangle loop, one should substitute in 
eq. ft4.8| ) the perturbative expression for the twist-2 part of the two-point correlator 

Ilnte 2 ) - nfV) = -L f 1 dyy n ln^f (4.11) 

57T Z JO ^ 

where \j? is the UV cut-off parameter for the composite operator 0{<9 Ml . . .d^ n }(p, and q 2 serves 
as a cut-off for the low-momentum region of integration. As expected, in the limit q 2 — > 0, one 
obtains a mass singularity. 

Substituting ( |4.11| ) into ( |4.8| ) and taking into account that {qi^ ■ ■ ■ qi^ n } differs from 
Qifix ■ ■ ■ Qifj, n only in inessential terms ~ Y 2 , it is straightforward to perform summation over n 
to get 

T B = r l dyln f^r^L eW -.„<„K> = /' dy ln (412) 



„o n 2 J tt 2 Y 2 Jo (p — yqi) 2 H 2 

Since (p—yqi) 2 = p 2 — 2y(pqi)+y 2 qf and (p— q\) 2 = —Q 2 , we can write (p—yqi) 2 = yp 2 + yq 2 — yQ 2 . 
Applying the SVZ-Borel transformation, we get 

^B { PT) = _ 1 f 1 Jy e -yQVyM> LmVM* ^ ( VV^\ 1 (413) 

2tt 2 Jq yM 2 \ \ /i 2 j j K J 

Comparing this result with the exact expression fl4.6|) , one can observe that the non-analytic terms 
of two expressions coincide. Hence, their difference, i.e., the coefficient function of the SD regime 
(see Fig.[7^, Figs.|8],|9] ) does not have non-analytic terms (mass singularities) in the q 2 — > limit 

(cf. 



For the "(0 2 ) -condensate" correction given by the diagram shown in Fig.[U]<i, the contribution 
corresponding to the lowest twist operators exactly reproduces the singular term ( f4.7| ) ~ l/(? 2 - 
Thus, this singularity will not appear in the modified OPE suitable for the nonsymmetric kine- 
matics q 2 Q 2 . 



23 



As we stressed above, the two-point correlator in (|4.8| ) is responsible for the long-distance 
contribution ~ and is not directly calculable in perturbation theory. However, we can write 

it through a dispersion relation 



1 r°° 5U (s) 

UJql) = - / — ^fds + ( subtractions ), (4.14) 
TV Jo s — qf 



where 5U n (s) = (Jl n (s + iO) — H n (s — i0))/2i is the relevant spectral density. As usual, we 
will model it by a "first resonance plus continuum" ansatz. A similar dispersion relation can 
be written for the perturbative correlator Il^ T (qf) ( |4.11| ), with the perturbative spectral density 
5il^ T (s) substituting the exact one. It is easy to realize that the ambiguity in the value of the 
/x 2 -parameter in eqs. ( (4.11|) , (|4.13|) corresponds to the ambiguity in the UV subtraction procedure 
for the correlators. However, the large-s behavior of the exact and perturbative spectral densities 
is the same, and it is not necessary to specify the subtraction procedure because the correlators 
appear in the modified OPE only through the difference U n (qf) — Ii^ T {q\) (see Fig.^|), which is 
UV finite. 

Writing explicitly the sum over the "hadronic " states in ( |4.14j ), we get 



» +1 n„( 9l ,y) = 'i" +1 n, w ..^,( 9l )y«^^'"- 



= fJOWfl + i r ,«£MW + ( subtractlons , (415) 

By analogy with the definition of the matrix element for the 7r-meson, we define {0\j(0)\4>,p) = if^. 
We can introduce the twist-2 distribution amplitude (f)(y) by treating the matrix elements of 
composite operators as its moments: 

<0,p|0(o)(r«9)"0(o)|o> = i n (Y P r(-iU) !\My)dy- (4.16) 

Jo 

We also used a convenient shorthand notation (y n ) = Jq 1 (p(y)y n dy. The contribution of the higher 
excited states (continuum) in eq.( [4.16[ ) is approximated, as usual, by the perturbative spectral 



density 5H^ T (s) = — (I/871") Jo y n dy starting from the continuum threshold s^. 

The parameters of the model spectral density in eq.( |4.15|) , namely, the mass m^, the residue of 
the first pole fy, the moments of the distribution amplitude (y n ) and the threshold should be 
extracted from the auxiliary sum rule for the moments of the scalar meson distribution amplitude. 
In this case, the additional terms in the modified OPE will decrease with increasing q 2 . As a result, 
the modified SR will reproduce the original SR for the three-point correlator valid in the region 
where both q 2 and Q 2 are large (see Fig.|T0D. Treating these parameters as known, we substitute 
( |4.15| ) in ( |4.8| ) and define the "bilocal" contribution in the r.h.s. of the SR for the three-point 



correlator. As a result, all the additional terms in the OPE (see the second row of Fig.|5]) can be 
written in the following form: 
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Now we can write our final expression for the r.h.s. of the modified sum rule: 
$(g 2 , Q 2 , M 2 ) = $ (a) + $ (6) + $ (c) + + A$ 



1 1 



o yM 2 



mi + ^ 



2tt 2 



■In 



(ScP+q 2 )y q 2y/M 2 



M 2 



+ 



2 ) 



Q 2 M 2 



(4.17) 



(4.18) 



It is straightforward to observe that this expression is well-defined for q 2 = 0. Below, in our analysis 
of the QCD sum rule for the F 7 * 7 .^ 7r o(g 2 , Q 2 ) form factor we will follow a similar strategy. Of 
course, instead of the toy scalar meson, the p-meson will play the role of the lowest state in the 
bilocal contributions. 



5 Mass singularities in the QCD case 



As discussed above, the B-regime for the correlator ( |2.^| ) corresponds to a situation when only the 
points Y and are separated by a small interval (see Fig.||]&). In simple cases, the short-distance 
coefficient function is given by a single quark propagator or a product of propagators. The long- 
distance contribution is represented by a particular two-point correlator of the electromagnetic 
current J M (A) and some composite operator (see E 
two lowest twists appear in this expansion. 



As we will see, only the operators of 



5.1 Terms with single-propagator coefficient function 

The simplest case is when the coefficient function is given by a single quark propagator S(Y) 
Y /2n 2 {Y 2 — iO) 2 . The relevant contribution (Fig.|6]5) can be written as 



2vr 2 ^ -n 



OO 1 



x 



S v0our j d 4 Xe~^ x (0\T{J^X)u(0)(d^...d^halM0)}\0) 

+ieu/3aaj d 4 Xe~^ x (0\T{J^X)u(0)(d^...d^)7 ff u(0W) 
+ { S v(3aa f d 4 Xe-^ x (0|T{J M (X)«(0)(^ 1 ...a Mn )T<x75«(0)}|0) 
+ie v(3aa f d A Xe~^ x (0|r{J M (X)«(0)(9 M1 ...9 Mn ) Tff «(0)}|0> 



(5.1) 
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where S vj3aa = {g v pg aa - SW 9p<r + 9ua g a p)- 

Let us consider the bilocal correlators with the right-sided derivatives 

R 5 n (q h Y) = J d*Xe-^ x (0|T{^(X)^(0)(F9) n 7 CT 75n(0)}|0), 

Rn{qi,Y) = J d 4 Xe~^ x (0\T{J^X)u(0)(Ydr 1(7 u(0)}\0). (5.2) 

The bilocals L\(q\,Y) and L n (qi,Y) with the left-sided derivatives can be treated in the same 
way. 

For any n, we can expand the current with derivatives over a set of traceless operators. More 
precisely, one should deal with traceless combinations of the indices /3, /ii, . . . , \i n . Therefore 
the essential contribution to eq. (|5.1|) comes from two types of operators, viz., the lowest twist 
operators which correspond to the traceless ■ ■ ■ , /i n } combination and the next-to-leading 

twist operators which contain one contraction of ~ or ~ g^. type. The operators with higher 
twists are accompanied by the factors (Y 2 ) 2 , (Y 2 ) 3 , etc. which cancel the singularity of the quark 
propagator 1/Y 4 and, hence, do not produce mass singularities. 



5.2 Factorization of the perturbative term 



The factorization procedure for the perturbative loop is illustrated on Fig.|7|a, b, c. The diagram 
0c corresponds to the expression ( |5.1| ) with the two-point correlators given by their lowest-order 
perturbative form: 



Rn(.qi,Y) 



12 



-iYk) n [k a kp - k a q ip ] 
k 2 {k — qi) 2 



(5.3) 



where d D k = d D k/ (2ir) D , D = 4 — 2e; here and in the following we use dimensional regularization 
for the UV divergencies and the MS subtraction scheme. 

To extract the contribution of the lowest two twists it is sufficient to keep only the terms up to 
Y 2 in the expansion of the integral ( |5.3|) in powers of Y 2 . Indeed, the twist- 2 contribution in ( |5 . 1| ) 
is obtained by taking formally Y 2 = in the numerator of the integrand. Terms proportional 
to Y 2 give the contribution of the twist-4 operators. Using this expansion (see Appendix B) and 
( |5.1| ), it is possible to perform summation over n and to integrate over d 4 Y. After simple but 
lengthy calculations we find: 



$ B(PT \q 2 ,Q 2 ,M 2 
x/e^/^ ln 



7T JO 



' Q 2 v/M 2 fj 1 



dye 



<fyy 



1 4 2 

2 <ry qjr 

M 2 M 4 



2M 2 

+ 2^ 

M 2 



4„,2 



3q y 



2 M 4 



(5.4) 



As expected, the terms proportional to q 2 In q 2 (given by the twist-2 operators) and q 4 In q 2 (pro- 
duced by the twist-4 operators), coincide with the non-analytic terms in 
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5.3 Factorization for terms proportional to the gluonic condensate 

In a similar way, we can consider the diagrams proportional to the gluon condensate (see Fig.|3|). 
In this case, it is convenient to perform factorization diagram by diagram, because different groups 
of diagrams correspond to different coefficient functions (CF) of the B regime. Applying the same 
technique that was used for the perturbative contribution, we get the following representation for 
Fig#,c: 



111 



(5.5) 



9M 6X tt 'Jo \q 2 y \y y 2 1 M 2 



i"7JSM —;*-)■ (5 - 6) 



The factorization procedure for these diagrams can be performed just like it was done in the 
perturbative case. The only difference is that now we should take the {GO) part of the correlator 
( |5.2j ) rather than its perturbative part. For Fig.p?, the bilocal contribution is 



I yq 2 [i 2 6 



and for Fig.^c we obtain 

* B M\Q\M 2 ) = ^ GG )^- 2 (5.8) 

18 x tt ' M 6 Jo y \yq 2 \y y 2 ) J 

Note, that the terms proportional to 1/q 2 are due to the traceless combination of indices 
/3,/ii, • • • , fi n in (|5.1j ), whereas the terms proportional to lng 2 correspond to the next-to-leading 
twist operators. 

Diagrams g can be treated in a similar way, their total contribution being 

* f+9 { q 2 ,Q 2 ,M>) = -y^GG)^ J^ye-^ (5.9) 
J (2y 2 + y)f yq y M , y^_^(yq 2 \ n ^{n + l)\ (1 — 2y) M 2 \ 



y 2 V M 2 ^ \M 2 I n\ J y 2 q 2 
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In this case, the short-distance part of the relevant bilocal contribution is given by a product of 
two quark propagators (see eq. (|6.36|) below). Extracting the (GG) contribution, we get 



*f +g (M,M 2 ) = -^GG)^J Q dye-*"" 



1 /W 2 y 2 g 2 J 

where (...) denote terms regular for q 2 = 0. 

5.4 Factorization for quark condensate corrections 

The factorization procedure should be applied also to the diagrams with the quark condensate. 
The relevant contributions are analogous to the power corrections of the toy scalar model. 

Consider first the diagrams with a soft gluon (Fig.f|). The first three diagrams (Figs.ffla, b, c) 
produce a p 2 -independent contribution into F(p 2 , qf, g|) which vanishes after the SVZ-Borel trans- 
formation in p 2 . The last three diagrams (Figs.§g, h, i) are regular in the q 2 — > limit, and there 
is no need for an additional factorization. Finally, the contributions of the remaining three dia- 
grams (Figs.f|(i, e, /) has a manifestly factorized form (see (|5.1|) )( jQ6|) ), i.e., they are completely 
absorbed by the bilocal terms in the modified OPE. 

For the diagrams with a hard gluon exchange (Fig.[|), the situation is very similar. The first 
three diagrams (Figs.[5|a, b, c ) produce a manifestly factorized contributions into F(p 2 , qf, q 2 ). The 
basic difference between them is that the coefficient function of the I?-regime is given by a product 
of two quark propagators while those related to the diagrams |5]&, c is formed by a product of two 
quark propagators and a gluon one. The contribution of Fig.[5|d is regular in the q 2 — > limit, and 
there is no need for an additional factorization. 

As mentioned earlier, the remaining diagrams |5|e-r do not contribute to F(p 2 , qf, q%) in the 
large-g 2 regime. However, a more careful analysis (see Section 7) of the bilocal objects corre- 
sponding to some of these diagrams shows that there may appear nontrivial contributions into 
F(p 2 ,qf,q 2 ) due to the so-called contact terms JH]. Within the QCD sum rule method, the 
contact terms, in particular, play a crucial role in establishing the correct normalization for the 
electromagnetic form factors at zero momentum transfer E3, 12(3 . 



After application of the subtraction procedure illustrated on Fig.|^, all the infrared singularities 
of the original sum rule cancel with the corresponding singular contributions from the diagrams in 
which the B-term is extracted explicitly. In our case, it is sufficient to consider the bilocal objects 
corresponding to operators of the lowest two twists. The terms remaining after this subtraction 
are regular in the q 2 — > limit and contribute to the modified sum rule. 



6 Bilocal correlators 
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6.1 Bilocals related to a single-propagator coefficient function 

Let us start with the case when the coefficient function of the 5-regime is given by a single quark 
propagator. Its contribution into the three-point amplitude T allv {(i\, p) can be written as 

r Y^d^Y r 

K^P) ~ J eVPY jY^W J e ^ lX ( l r ( J ^ X )^(°)^^7^(V))|0)rf 4 X. (6.1) 

As emphasized before, the correlators (bilocals) (0\T{J fM (X)O^\0)}\0) of the J M (X)-current with 
composite operators 0^(0) are sensitive to the long-distance (~ dynamics and, for this 

reason, they are not directly calculable in perturbation theory. The standard way out is to write 
them in a dispersion form and assume the simplest ansatz ("lowest resonance + continuum") 
for the model spectral density, with the continuum starting at some effective threshold s p . The 
parameters of the model spectrum can be extracted from an auxiliary QCD sum rule. The 
quantum numbers of the electromagnetic current J M (X), which appears in all the bilocals, dictates 
that the lowest resonance should be represented by the p°-meson []. In particular, for the bilocals 
corresponding to the-single propagator coefficient function, we can write 

„ . T ,. 1 f°° 5R n (s) , . . 

R n (qi,Y) = — / T ds + ( subtractions ) 

7T Jo s — qf 

1 f°° 6R PT (s, Y) 

= "p°-meson contribution" H — / ds — \ h ( subtractions ), (6.2) 

7T Js p s — qf 

where SR n (s, Y) is the relevant discontinuity: 5R n (s, Y) = (R n (s + iO, Y) — R n (s — iO, Y))/2i. 
Picking out the p°-meson term 

5 (+ V - ml) ^ E \ P ° X ; P) ( P ° x ; P \ c t, (6.3) 

in the sum over physical hadronic states (with A being the helicity of the p°), we extract the p° 
contribution. As a result, we obtain a set of matrix elements which can be parameterized as 

(0$(OM(Y)\P° X=0 -,P) =ip a fj4> v p {Yp, l i 2 ) + -.. (6.4) 

(0|^(0) 7 ^OniP|A|=i;P> =4 f p 7 m p <f> v p± (Yp, t i 2 )+ia vl p a tfm^Y) ^(Yp, p 2 ) + ■ ■ ■ (6.5) 

(01^(0)7^(^)1^1=!;?) =e aoi p 5 eip p Y s ftcj>*(Yp,p 2 ) + .-- . (6.6) 

Here only the twist 2 terms are written explicitly, and the dots stand for the higher twist contribu- 
tions, e a is the polarization vector of the p°-meson, and the helicity components have an evident 
interpretation in terms of the longitudinal and transverse polarizations: p x =o = Pl 5 P°\\=i = P±- 
In a standard way, the functions <p p (Yp, p 2 ) can be related to the usual wave functions (p p (y, p 2 ) 
describing the light-cone momentum distribution inside the p: 

$\Yp,p 2 ) = ^dye-^tpf&p 2 ), (6.7) 

4 Strictly speaking, we deal with an isotopic mixture of p°- and ^-contributions which degenerate in the chiral 
limit m u = md — to an effective state with the mass of the p-meson (cf. Q ) . 
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with /x 2 being the renormalization parameter for the relevant composite operators. The constant 
fixing the normalization of the simplest wave function is known from previous QCD sum rule 



studies: fY ~ 200 MeV |18], |48|], while the constants and a V i in eqs. (|6.5|), fl£T6| ) can be fixed by 
equations of motion (cf. E3|, |S4J), which form an infinite set of relations connecting the moments 
of different wave functions (see Appendix C). 

For our purposes, it is more convenient to write down the matrix elements in a form suitable 
for an arbitrary polarization of the p°-meson: 



+ia vl p CT f v m p (e^Y) [^(Yp^ 2 ) + C [vl] ,Y 2 ^ {Y p, u 2 ) + ■■■ 
+f V m pC[V2]4 (YM X) Y) - e^Y 2 /A) <f>™(Yp, /i 2 ) + • • 



(6. 



(0m0h^(Y)\Px,P) =e^ P e^PpY p [<p A p (Yp^ 2 ) + C A4 Y 2 <pf(Yp, f, 2 ) + ■ 

Since the C-parity of the p°-meson is negative, its wave functions have the following 
(here and below y = 1 — y): 



. (6.9) 
properties 



Jo 



dy Lp 



V,Vi,[V2]A,A,A4: 
P± 



(y) 



i 



dyy^ nvi] \y) = i- 



In the relations above, we have explicitly displayed wave functions up to twist 4. 
for a longitudinally polarized p°-meson, 



(6.10) 
Note that, 



_A=0 



iPa/m p + O (m p /p z 



as p z — > oo, and the twist-2 part in eq. (|6.8| ) coincides with the well known definition 

Let us consider first the twist-2 p°-meson contribution in eq. (|6.2j ). Applying 
(|6.9| ), we obtain: 



' OO 1 N 

Y.-M^y) 



-i){^m pl 



\n=0 



nl 



p° 



m 2 -q\ 



„_„ , dye-^y^Ky) 



-ia viqia e p {e*Y) I dy e^** <pY}(y) 



~M and 



(6.11) 



\n=0 



nl 



m — q 1 jo 



(6.12) 



Here we use the shorthand notation: 



A=0,±1 



-9fMT + 



mi 



(6.13) 
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Substituting eqs. ( |6.11| ), ( |6.12| ) into eq. ( p.l| ) and extracting the proper tensor structure we 



get: 



where 



4vr />, /•! , 1 



3 m \ — q{ Jo p 



A ,A 



(6.14) 



72 + yqif = -q\ yy + qlv + p 2 y 



is the virtuality of the hard quark written in the "parton" form. This formula includes an extra 
factor of 2 which appears when one adds the contribution of the correlators L n , L\ and uses the 
symmetry properties of the wave functions ( |6.10D . It should be noted that the wave functions <p^ ± , 
(p^* and y?j^ 2 ' 4 do not contribute to the form factor F which is considered here. 

Deriving, in a similar way, the twist-4 contribution and applying the SVZ-Borel transformation 
to the resulting amplitude, we obtain: 



4vr f v p m p 



3 m\ + q 2 Jo y 2 M A 



(6.15) 



x 



«vi f?m p {^(y) - ^ ) - (1 + 2y) U (y) - <%) 



yM 2 



yM 2 



As expected, the twist-4 contribution is suppressed by one power of 1/M 2 . 

Now, taking q 2 = and introducing the integration variable s = yQ 2 /y, we can write this term 

as 



A{Q 2 ) 



47r fj f 00 ds 



3Q 2 m p Jo M 4 



s/M 2 



rV I VI 

a V i f p m p \ <p p± 
s 



■ S + Q 2 , 



AC [vl]A {s + Q 2 ) [yi]4 



s + Q 2 



+ 



4C A4 (s + Q 2 j A4 



M 2 Q 2 
s \ 



M 2 Q 2 



s + Q 2 



,s + Q 2 
(6 



Note, that this representation has "wrong" powers of the SVZ-Borel parameter M 2 compared 
to the canonical form 



$(M 2 ,Q 2> 



ttM 2 jo 



p(s, Q 2 ) ds. 



Using the transformation 



ds e 



s/M 2 9{s)_ 

M 2 



ds e- s/M '' 



9M 
M 4 



ds e~ s / M \g(0)6(s) + g'(s)), 



ds e - s / M2 (g'(0)5(s)+g(0)5'(s)+g"(s)) 



(6.17) 

(6.18) 
(6.19) 



etc., we can always cast the terms with the "wrong" powers into the canonical form ( |6.17 ). 
In particular, taking the asymptotic forms (see Appendix C) 



a vi 



- f A 

40' Jp 



> f) P 

, vvi = <pvi = 6 °yy O - *)> <pa = <pa = 6 yy> 



(6.20) 
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as the simplest estimate for the twist-2 contributions, we get 

P P 2(s,Q 2 ) = 9L(s) = SvrUT iQ4 - 2 " Q P (6.22) 



Note, that the s-integral over the latter spectral density is zero. 

The asymptotic forms for the twist-4 distribution amplitudes can be directly extracted from 
the corresponding correlators: 

<P\yi]4 = <Pyi]A = 420(^) 2 (2y - 1), VA 4 = <p% = 30y 2 f- (6.23) 

From the equations of motion (see (C.4), Appendix C) it follows that C[vi]4 = 5Ca^/7- Taking 
into account this expression and eqs. (|6.23| ), we obtain: 



where 

g AiP (s) = -4tt 2 (/;) 2 40 C A4 { f + Q ^ . (6.25) 
Hence, the relevant spectral density is 



p»(s,Q 2 ) = -^{f V P f 40 <W (3g2 u + Q4) = 9'IM (6-26) 



\2\ _ Q2(rV\2 An ri r\2 u ^ ~ r ^ > " 



6.2 Continuum contribution 

As noted earlier, in the basic OPE for the small-g 2 kinematics, one always deals with the difference 
between an "exact" bilocal correlator R and its perturbative analog R PT (see Fig||). Since the 
ultraviolet behaviour of these two correlators is the same, there is no need to explicitly specify a 
subtraction prescription for the correlators. 

Now, incorporating our model for the bilocal correlators, in which the contribution due to 
higher excited states is approximated by the perturbative spectral density (see (|6.2j)), i.e., by the 
continuum starting at s p , we can easily write down an expression for the difference between the 
continuum contribution to R and the perturbative bilocal R PT . Then, substituting the result into 
the original expansion (|5.1| ) and performing some straightforward calculations, we obtain: 

-n/M^y e q 2 y/M 2 



^B{cont) _ ^B(PT) = 1 f 1 dy— e - Q2 
n Jo M 2 

( . « _|_ a 2 M 2 \ 

(6.27) 



W 2 ^ 8 ") + w[ q qSp + j) 

The terms collected in the ( ) brackets correspond to contributions due to operators with 
twist 2 and 4. Note, that these terms exactly cancel the logarithmic contributions q 2 In q 2 , q 4 \nq 2 
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present in the coeficient function of the unit operator for the usual OPE valid in the large-g 2 
kinematics. As a result, the non-analytic terms are replaced by the combinations q 2 In (s p + q 2 ) 
and q A In (s p + q 2 ), which are "safe" in the q 2 — > limit. On the other hand, for large q 2 , the 
usual OPE without additional terms must work, i.e., the difference between "exact" bilocal term 
and its perturbative analogue must vanish faster than any power of 1/q 2 in the large-g 2 limit. It 
is easy to check that the terms on the r.h.s. of eq. (|6.27|) behave like 1/q 2 when q 2 — > oo. To get 
the total expression for the additional terms, we should add the p-contribution to eq.( |6.27f) . The 
p-contribution also has the l/g 2 -behaviour for large q 2 . However, to produce a perfect transition to 
the pure SD-case, the l/q ,2 -terms must cancel. Basically, this means that using a rough model for 
the correlator one should not rely too heavily on the extrapolation of our result ( 6.27| ) beyond the 



region q 2 <m 2 . However, we can require that the model, at least, should provide the cancellation 
of the 1/q 2 terms. If we choose the asymptotic form for the lowest-twist distribution amplitudes, 
the cancellation of the l/g 2 -terms produces the estimate: 

s 2 p = 8AfJ) 2 m 2 p . (6.28) 



Our result ( p. 27\ ) simplifies in the q 2 — > limit: 

foB(cont) _ (f)B(PT) _ 1 f°° J -s/M 2 f_ 
~ 7TM2 Jo Q2 



Sp + 



M 2 p M 4 2 



(6.29) 



where s again is defined by s = yQ 2 /y. The twist-2 contribution can be represented in the 
canonical form with the following spectral density: 



( Q 4 - 2sQ 2 ' 
(s + Q 2 Y 



pi(s,Q 2 ) = -2s p ^^ ■ (6.30) 



For the twist-4 contribution in eq. ( |6.29|) , we obtain: 



Li i ,\2s 2n,2$ s2 6sQ 2 + Q 4 j a . . 

Pi {S,Q) = S p Q (g + Q2)6 = f4,cont(s), (6.31) 



where 



s 2 s 2 Q 2 

94,cont(s) = — (6.32) 



2 (s + Q 2 ) 4 

Using ( |6.28 ), the p-meson contribution from (|6.22|) can be rewritten as 



a. m= U.) = l®^%P-. (6-33) 

For the next-to-leading-twist contributions, similar duality arguments give the following esti- 
mates for the normalization constants of the relevant two-body distribution amplitudes: 

C7 A4 = — ^ 2.28 • 1(T 2 GeV 2 ; C wm = — 1.63 ■ 10" 2 GeV 2 . (6.34) 
60 84 
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Hence, the expression (|6.26|) can be written as 



,S(, Q ^-2^ (3S2 : 6 l: g4) (6.35) 



3my (s + Q 

If the distribution amplitudes deviate from their asymptotic forms, then the duality condition 
between the p-meson and the continuum should keep the form of the expressions ( p. 33] ), ( |6.35| ), 
but instead of s p we should get an effective duality interval. 

Note that for fj = 0.2 GeV, our estimates give s p ~ 1.37 GeV 2 for the duality interval. This 
value is in good agreement with the "canonical" one: s L p D = 47r 2 (/^) 2 ~ 1.58 GeV 2 . The latter 
follows from the local duality considerations for the two-point correlator of two vector currents. 

6.3 Twist-3 bilocals for two-propagator coefficient functions 

Next in complexity is the contribution related to the coefficient function formed by a product of 
two propagators S(Y, Z) and S(Z, 0): 

F?$ = — / d A Y e^ Y d A Z Z ~ Z) \, — t~~t E -T-T Yfl1 ■ • • Y " n ZVX --- ZVm ( 6 - 36 ) 
3 J 2tt 2 (F-Z) 4 2vr 2 Z 4 n j^ =0 n\m\ V ; 



x 



J d A Xe~^ x (0|r{J M (X)u(0)(d m ...^ • 



Here we explicitly extracted the bilocal correlator containing a composite operator composed of 
two quark and one gluonic field. Note that the gluonic field A b (Z) here may be treated as taken 
in the Fock-Schwinger gauge, i.e., it can be substituted by 

A\{Z) =Z v J 1 a G b ^(aZ)da . (6.37) 

As a result, the p°- meson contribution is determined by the following matrix elements: 

(0\u(Z 1 ) lpl5 g s Gl y (Z 3 )MZ2M;P) = Vpe^enPee^ f£ p ^ p {Z lP) p 2 ) 

+higher twist contributions (6.38) 

(OluiZ^igsG^Z^MZ^P) = p p (j^eW-^eW) fj?) 

+higher twist contributions . (6.39) 

In a standard way, we can introduce the momentum distribution amplitudes^ <f3 p A (yi)'- 

4>l A (Z iP ,^) = f\dy] 3 ^f'S*^). (6.40) 
They have the following symmetry properties: 

(P3 P {Vi, V2] ya) = vtfa, yr, y 3 ) , (fspivu Vs) = - fX p {V2, Vi\ ys)- (6.41) 

5 [dy} 3 = dyidy 2 dy 3 S(l - J2iVi) 



34 



In our definition, the normalization constants f 3 = 0.6 • 1(T 2 GeV 2 , f 3 v = 0.25 • 10~ 2 GeV 2 [4 
are factored out, so that the distribution amplitudes are normalized to unity: 



Following the procedure described in Sec. 2.1, we find the p°-meson contribution: 



$ B (2) 



?7T f V p m p rl rl rl 



3 m? + q 2 Jo 



where 



da a dp [dy] 3 e b/aM 
Jo Jo 



x {fspfipiyiiyxys) 



Cl 



d\ 



a 2 M A 2a 3 M 6 



C-2 



do 



a 2 M 4 2a 3 M 6 



(6.42) 



(6.43) 



a = a(3y 3 + y 2 , 



q (a (3y 3 + 2a(3y 2 y 3 - ot(5y 3 + y 2 - y 2 ) + Q (a(3y 3 + y 2 -l) 



and 



Cl 

d\ 

c 2 
d 2 



(6.44) 



(a(3y 3 )/(a(3y 3 + y 2 ), 

ci ( ~q 2 (a 2 (3yl + 2ay 2 y 3 + a(3y 3 + y\ + y 2 ) + Q 2 
(a(3y 3 + 2(3y 2 )/ (a(3y 3 + y 2 ), 

ci ( q 2 {a 2 (3y 2 3 + 2af3y 2 y 3 + af3y 3 + 2/?2/| - y 2 + y 2 ) + Q 2 (1 - 2(3) ) . (6.45) 



The perturbative spectral density for all of these correlators is suppressed by 0(a s /7r)-factor, 
and for this reason we neglect here the contribution due to higher states. 
Taking q 2 = 0, we get 



and 



a = a(3y 3 + y 2 , b = -Q 2 (1 - a) 

ci = (a(3y 3 )/a, d x = c x Q 2 \ 

c 2 = {af3y 3 + 2f3y 2 )/a, d 2 = d Q 2 (1 - 2(3). 



Introducing the variables: 



n2n w n2 1 - 2/2 - aPV3 

s = Q (1 - a)/a = Q ■ — , 

y 2 + a(Jy 3 

and integrating over (3, a, we obtain the representation 



v(s) 



Q 2 



s + Q 2 



(6.46) 



(6.47) 



(6.48) 



p 



87T f, 



v 



3 m p Jo 



ds e~ s ' M2 



{s + Q 2 ) [K, , a( , , ^ (v(s) - y 2 ) 

W" Jo [dyh 9im - V(S)) ~~d X 
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fs P fipiVh V%\ Vz) (V2 + V3- v(s)) 



[s + Q 2 
2M 2 



(6.49) 



(s + Q 2 ' < 



V2 + y?,-v(s) -2y 2 hi 



2M 2 



V2 + V3- v(s) + 2(v(s) - y 2 ) In 



'v(s) - y 2 
'v(s) -y 2 

. V3 , 



To estimate these contributions, we use the asymptotic forms of the corresponding three-body 
p-meson distribution amplitudes [@: 



WAivmys) -> (PaAivmys) = 360y!y 2 yl, 
(Psvivmva) -» VsMfiStel/s) = 7! (1/1 - y 2 )y\y 2 y\. 



(6.50) 



6.4 Twist-2 bilocals for three-propagator coefficient functions 

The bilocals associated with the coefficient functions given by a product of three propagators 

— 2 

can appear in the (ipip) quark condensate diagrams of the unmodified OPE (see Figs.^a — r). 
Furthermore, for large and moderate qf, only the diagrams [5|a — d contribute to the invariant 
amplitude F we are interested in. So, let us consider them first. In fact, among these diagrams, 
only ^ and §c produce bilocals with the three-propagator coefficient function. After some algebra, 
we obtain 

m _ 32n 2 a s (uu) r . Y ( pA Y p « 1_ 
^ ~ 27 J \p 2 J 8n 2 Y 2 n\ ' 

x f d 4 Xe~^ x (0|T{J M (X)n(0)(^ 1 ...^J 7 ,7/375«(0)}|0). (6.51) 



Adding the charge conjugate contribution produces an extra factor of 2. The correlator which 
appeared in eq. ( |6.51| ) can also be treated as a distribution amplitude <p 7 * (y, <? 2 ) of a photon with 
virtuality q\ = —q 2 . The perturbative spectral density for this correlator is zero, so the natural 
approximation is to model <p 7 * (y, q 2 ) by explicit contributions from lowest resonances: 



M^ 2 ) = ^4rl^M + *^T^* T M + ■■ ■ , (0.02) 
m P + Q TUpi + q 



where the p-contribution is determined by the matrix element 

(0m0)a u ^(Y)\p° x ;P) =i{e<f> p p - ef p v ) f p T p (Yp, p 2 ) + higher twists . (6.53) 

Here, as usual, a v p = f [ 7i ,, 7/3], p = qi + qi- For the SVZ-transform $7* (3) we thus get: 

r(3) _ 64tt 2 a s (qq} f 1 f 1 , , „ 0tf v g-Q* v )Mi-v0)M>) P <Pj* fa <?*) , fi ^ 

27 Af6 J0J0 V P (1 - ' 1 J 
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The distribution amplitude of the real photon (p~,(y) = y2 7 *(?/,g 2 = 0) was analyzed in ref.|27|]. It 
was argued there that <p 1 (y) has the "asymptotic shape": 

2k 



mi 



(6.55) 



with k ~ 1.3 and 7 (y) = 6y(l — y) (0 7 (y) is the normalized photon distribution amplitude i.e., 
its zeroth moment equals 1). Taking q 2 = and introducing the variable s = Q 2 yj3/(1 — y(3), we 
then get: 

d> 7(3) = ™? ka A m) l r e -^ 9l (s)ds. 



27 M 6 Q 4 m 2 Jo 



where 



s/s+Q 2 



y 



(6.56) 



(6.57) 



For the model </> 7 (y) = 6y(l — y) advocated in p7f , we obtain 



g 7 (s) = 6s In 



Q 2 + s Q 2 



Q 2 



(6.58) 



Using the formulas ( |6.18D ,( |5'.19p , we can convert the expression ( p.56| ) into the canonical form. 
Note, however, that the second derivative of g 7 (s) contains the 1/s-singularity, and one should 
be careful when calculating the relevant spectral density p 7 (s, Q 2 ). The simplest procedure is to 
represent the In s term as lim^^g ln(s + A 2 ). Then application of eq. (|6.19| ) is straightforward, and 
we get 



= i2g^ a ,m 2 lim 

27 M 2 Q 4 m 2 x^o 



1 ^ 9 

In- -2 



-s/M 2 



' s 2 + 3sQ 2 + AQ 4 



(s + Q 



2\3 



s + X 2 



ds 



(6.59) 



This trick amounts to using an alternative regularized form for the (l/s) + -distribution usually 
defined by 

s+ Jo s 



-ds. 



(6.60) 



7 Bilocals and contact terms 

A special care must be taken about the correlators containing the Dirac operator 7 M -D M acting 
on the quark field ip. Since the correlator is a T-product of the electromagnetic current and 
a composite operator, applying the equation of motion one gets the 5^ (X)-function, i.e., the 
external vertices of the bilocal are contracted into a single point and it reduces to a g 2 -independent 
constant. 



Let us sketch a simple derivation for such terms (see, e.g., |20|). Using the functional repre- 
sentation for the correlator 

(0|T{...^(X)W(0)}|0)= I V[i>]V[^]V[A]{..4{X)Vm}^v{i I C(Z)d 4 z), (7.1) 
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where C(Z) = ip(Z) iVip(Z) + • • • , we can write 

6 



W(0) exp [ij £(Z)d 4 Z^J = --— exp (i| £(Z)d 4 z). (7.2) 
Integrating by parts in ( |7.1|) results in 



(7.3) 

It is the derivative 5ijj(X) / 5ip(0) that produces the 5^(X) term mentioned above. 

The contact terms play an important role in all applications of the QCD sum rules to low- 
momentum behaviour of hadronic form factors. In particular, without them, it is impossible to 
satisfy the Ward identities fixing the pion form factor normalization at zero momentum transfer 



7.1 Separating short- and long-distance contributions 

Consider the hard gluon exchange diagrams shown in Figs.|5|a — r which produce, in the 5-regime, 
the bilocals associated with the three-propagator coefficient function. Take first the diagrams 
|]e,/. Their contributions are 

32z 32i 

~ ^2^,2^4 1l^ e oivqiq2 1 ^ J ~ ^4„2„2 lift € avqiq 2 i I'-^J 

P H\H2 P 

i.e., they do not contribute to the invariant form factor F. However, as we will see below, in 
the B-regime, equations of motion "extract" the appropriate tensor structure, i.e., these diagrams 
cannot be ignored. The relevant term from the three-point correlation function can be written as0 

T B ' 5ef (a, ao) = _ 3271-2 a -s(uu) r 4 ipZ 1 ~ 1_ 



n=0 



X 



J d 4 Xe-^ x (0|r{J M (X)«(0)(5 M1 ...^H0)}|0). (7.5) 



Extracting the bilocal term from ( |7.5| ), one should pick out the traceless combination of indices 

fj,i, . . .fi n , i.e., the lowest-twist term which gives the leading power contribution with respect to 
l/p 2 , l/qi- Introducing the notation 

IW.^CgO = J d 4 Xe~^ x (0\T{ J,(X) u(0){d^ . . . ^>(0)}|0>, (7.6) 

we can represent the correlator (|7.6|) in the following form: 

n A4{Ml ... Mn} (gi) = A^ n \ql) q lfl { qilll . . . g lMn } + {q lfl , . . . g lMn } + 

+C (n \q 2 1 )g,{, 1 q 1 , 2 ...q 1 , n }, (7.7) 



3 To get the total contribution, one should add also the charge conjugate term. 
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where {. . .} denotes the traceless-symmetric part of a tensor. Because of the electromagnetic 
current conservation, we have the constraint gi^ n M { Atl .., Atn }(gi) = which produces some relations 
between the invariant functions A^ n ' , B^ n \ C^ 1 ' . Using the formula from [f49| 



n + 1 
2n 



xQifn ■ ■ ■ Qi^A > 



we obtain 



' A (n) , (n + 2) " 
, + 2(n+l) , 



+ C (n) = 0. 



Contracting ( |7.6|) with gives 

n^^Cft)^ = (>g 2 ^i^ + c^ ( 



n + 1 



n 



(7.8) 
(7.9) 

(7.10) 



Furthermore, applying the technique symbolized by eqs. ( |7.1| ) - ( |7.3| ) (see Appendix E) we obtain, 
in the leading-twist approximation: 



\n— 1 



{9 



1/^2 



n M{Atl ... Mn }(gi)fl' wl ^ (-i 

-2<im) - \q le J d 4 Xe~^ x (0|T{ J M (X) «(0){9 W . . . ^ n }^(0)}|0) 



(7.11) 



where the first contribution inside the brackets is just the contact term, while the second one 
contains the correlator related to the photon distribution amplitude ip y *(y,q 2 ) ( |6.52| ): 

3 



Vl^l-'-Mn} 



{<llH2 ■ ■ ■ Ql^nf 



-2(uu)--q 2 I dyy n - L ip r (y,q 2 ) 



(7.12) 



Note that the term containing <^ 7 *(?/,g 2 ) vanishes for q 2 = 0. Using the formulas from Appendix 
D, one can easily perform the necessary contractions: 



+ 



+ 



n + 1 

C(n) 



-^r 71-1 C^-i(^) - qifj, — i~ n ~ 2 cl_ 2 (v) 



(7.13) 



where C^(rf) are the Gegenbauer polynomials and the notation n = i (qiZ) / \J —Z 2 q 2 , r = 

—i\J—Z 2 q 2 /2 is introduced. 

The tensor structure (~ p a €^ uqiq2 ) we are interested in, can be only produced in (|7.13f) by the 
terms ~ Z^. Hence, other terms can be ignored. Combining now eqs. (|7~9|) - (|7.12|) , we get, 
modulo the next-to-leading twist contributions: 



Z^-i q x Z) 



2 



n-l 



2{n + l) 
2n 2 



+ 



n 



(7.14) 



(n + l)(n + 2) 



-2{uu) - -q I dyy n <p r (y,q ) 
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Substituting ( |7.15| ) into (|7.5|), integrating over d 4 Z and summing over n by using the generating 
function technique we get for the SVZ-transform of the contact terms: 



$5e/(C) = 256tt 2 a s (qg) 2 
27 Q 2 M 6 



(7.15) 



1 daaipdy *"^ ~y® ifAvA 2 ) e^T^m-y^u* 



The non-contact terms give 

27 Q 2 M 6 4 y 

where y/3 = l—y/3, (3 = 1 — (3 . Analyzing the remaining bilocal contributions capable of producing 
a coefficient function of the three-propagator type (see diagrams of Fig.|5| b,c,i,l,o,p), we found 
that they do not contain the contact terms with the tensor structure p«e M ^ gig2 . Furthermore, 
there are no contact terms in the bilocals corresponding to one- and two-propagator coefficient 
functions. 

For q 2 = 0, the contribution of the non-contact terms vanishes while the contact terms give: 



$ C (Q 2 ,M 2 ) 



2567i 2 a s (qq) 2 r°° 
27Q 6 M 6 Jo 



-s/M 2 



ln!±£-2. Q2 



Q' 1 



sds. 



(7.17) 



Representing again Ins = lim A 2^ hi(' s + A 2 ) an d using eq. (|6.19| ), we obtain 



$ C (Q 2 ,M 2 



256n 2 a s (qq) 2 



lim 



27Q e M 2 \ 2 ^o 



In 



A 2 



3 + 



s/M 2 



's 2 + 3sQ 2 + 6Q 4 

, ( S +Q 2 r 



s + \ 2 



ds . 
(7.18) 



8 QCD sum rule in the small-g 2 kinematics 

Collecting now all the contributions, we obtain the theoretical part (the modified OPE) of the 
QCD sum rule for the form factor F 7 » 7 *^ 7r o(g 2 , Q 2 ) (see Fig||): 

$(g 2 , Q 2 , M 2 ) = $^(g 2 , Q 2 , M 2 ) + $£ G> + ^ + 



+ $ P(2) _ $ B ^ _ qB _ $ B + 

+ $7(3) _ $B _|_ r^5e/(C*) + <£5e/( 7 *) 



where the first two rows correspond to the original OPE valid for symmetric kinematics. Each of 
the next rows represents the additional terms corresponding to different types of the coefficient 
functions. As explained in Section 4, all the terms of the standard OPE, which are non-analytic in 
the q 2 — > limit, are cancelled by the corresponding B-contributions. As a result, the coefficient 
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functions of the SD- regime are analytic functions of q 2 (cf. [0). Substituting explicit expres- 
sions for all the terms which appear in eq. flS.ip gives the following expression for the SVZ-Borel 
transformed OPE for the three-point correlator valid in the region of small q 2 : 

$(g 2 ,g 2 ,M 2 ) = -^|/ o 1 ^ e-^y {(l + ^ggiW) + 



+e 



,g 2 y/M 2 



2y ( 2 , 



+ q 2 )y _ \ y 2 



M 2 



(s P + q 2 )y 2 , s 2 p 



M 2 



-Q s p 



9 7T 



2M 2 Q 2 M 4 Jo y ^ n! \M 2 



n-l 



+ 



64tt 3 , v 
+ -7T^a s {qq) 



243 



- v2 g 2 64vr 3 2 
—Us(qq) 



Q 4 M 2 



27 



2Q 2 M 4 



4tt 2 /> p /-i 1 
dy— 



x 



3 m 2 + q 2 Jo y 2 M 2 
4Ci 



yM 2 



yM 2 



x \f£p<p£ P (vi,y2',y3) 



3 m 2 + q 2 Jo 
Ci c?i 



a 2 M 2 2a 3 M 4 



data d(3 
Jo Jo 

■ fsp'PspiyuV^ys) 



b/aM 2 



C 2 



do 



a 2 M 2 2a 3 M 4 



" 27 M 4 dyd/? (l-y/3) 3 



256tt 3 a s (gg) 2 yyi 0$ - P)y p y0( q 2 fj-Q 2 )/«i-yp)M 2 ) 
27 Q 2 M 4 yy ^ (l-y/5) 3 



647r 3 a s (gg) 2 
■9 



(8.2) 



9 Q 2 M 4 7 77o ^ * (l-ya/3) 3 ^ 
This expression can be used as a starting point for constructing sum rules for the form fac- 
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tor F 7 . 7 . 7r o(g 2 , Q 2 ) at q 2 = 0, 



its derivative ■^F 1 * 1 *^{q 2 , Q 2 



or even for studying the q 2 



dependence of F 7 . 7 .„-o(g 2 , Q 2 ) in the region of small q 2 <m 2 for fixed values of the large virtuality 
Q 2 >2GeV 2 . 

For the constants present in eq.(|8.2j), we use the following numerical values: f] 



0.2 GeV, m p = 0.77 GeV ; the constants f p = —f^rrip/A, ayi = 1/40 are obtained from the 
equations of motion (see Appendix C), the values f£ = 0.6 • 10" 2 GeV 2 , f¥ p = 0.25 ■ 10~ 2 GeV 2 
are taken from the QCD sum rule estimates given in ref. [fig] , and the value k = 1.3 for the photon 
distribution amplitude ( |6.52|) was taken from ref. [27] . The quark and gluon condensate values are 
standard: ({a s /7i)GG) = 0.012 GeV 4 , a s (qq) 2 = 1.8 • 10~ 4 GeV 6 . For the continuum threshold in 
the p-channel we take the standard value s p ~ 1.5 GeV 2 obtained from the QCD sum rule for the 
p-decay constant fY [pl|. To estimate relative importance of various contributions, we take the 
asymptotic forms for y? 3 y(2/i, 2/2, 2/3), ^3.4(2/1,2/2,2/3) (eq.( p.50| )) and <pvi(y), (Pa(v) (eqs.(|C77|),(P^)). 

Our numerical analysis shows that the most important contributions in ( |B.2| come from: a) 
SD-regime (first five rows of (|8.2j)), b) term with the nondiagonal correlator (photon distribution 



amplitude (Pj(y)) and c) p°-meson contribution with the leading twist wave functions (B-regime) 
in the diagonal correlator. The terms associated with three-particle twist-3 wave functions are 
small: their contribution into the sum rule is of the order of a few percent. We evaluated the 
terms corresponding to the next-to-leading two-particle wave functions (of twist-4) and observed 
that they are also suppressed. The contact-type power corrections are small as well. 

Below, we consider only the simplest sum rule for F 7 » 7 . 7r o(g 2 = 0,Q 2 ), keeping in it only the 
most important terms listed above and contact terms. All the necessary expressions substituting 
the terms from eq.(|3.2|) by their q 2 = limit were given in the preceding sections. Combining 



them together, we obtain the QCD sum rule for the 77* 
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form factor: 



Q 4 -6sQ 2 + 3s 2 (s 2 



7T Ol 

+ — (— GG) 

9 7T 




1.3) 



It is understood that the sum rule is taken in the limit A 2 — > for the parameter A 2 specify- 
ing the regularization which we used to calculate the integrals with the (l/s)+ distribution in 
eqs.( 6.59f) , (|7.18 ). Furthermore, in this sum rule, we model the continuum by an effective spectral 
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density p e ^(s,Q 2 ) rather than by p PT (s,Q 2 ), including into p e ^(s,Q 2 ) all the spectral densities 
which are nonzero for s > 0, i.e., pf(s,Q 2 ) flOPP , p p 2 (s,Q 2 ) flQgp , pf(s,Q 2 ) flOlQ , p£(s,Q 2 ) 



( |6.35| ), p 7 (s, Q 2 ) ( |6.59| ), p c (s,Q 2 ) ( [7.1 8| ) and also an analogous contribution from the gluon con- 



densate term. 

We studied the stability of our sum rule with respect to variations of the SVZ-Borel parameter 
M 2 in the region M 2 > 0.6 GeV 2 . A good stability was observed not only for the "canonical" value 
Sq pa 0.7 GeV 2 , but also for smaller values of Sq, even as small as OAGeV 2 . Since our results are 
sensitive to the so-value, we incorporated a more detailed model for the spectral density, treating 
the yli-meson as a separate resonance at s = 1.6 GeV 2 , with the continuum starting at some larger 
value sa- The results obtained in this way are very stable with respect to variations of the SVZ- 
Borel parameter M 2 and, for M 2 < 1.2 GeV 2 , show no significant dependence on sa- Numerically, 
they practically coincide with the results obtained from the sum rule ( |8.3j ) for sq = 0.7 GeV 2 . 

In Fig.|TT], we present a curve for Q 2 F 77 . 7r o(Q 2 )/47r/ 7r calculated from eq. (|8~3]) for s = 0.7 GeV 2 
and M 2 = 0.8 GeV 2 . One can see that it is rather close to the curve corresponding to the Brodsky- 
Lepage interpolation formula ^/^^^.^(Q 2 ) = l/(l+Q 2 /47r 2 / 2 ). It is also close to the curve based 
on the p-pole approximation tt f n F(Q 2 ) = 1/(1 + Q 2 /m 2 ). It should be noted, however, that the 
closeness of our results to the p-pole behaviour in the Q 2 -channel has nothing to do with the 
explicit use of the p-contributions in our models for the correlators in the g 2 -channel. Recall, that 
we take q 2 = there, hence, these correlators simply specify some constants. The basic reason for 
the Q 2 -dependence of the p-pole type is the fact that the pion duality interval s ~ 0.7 GeV 2 is 
numerically close to m 2 ~ 0.6 GeV 2 . 

For Q 2 < 3 GeV 2 , our curve goes slightly above those based on the p-pole dominance and 
BL- interpolation (which are close to the data ||). This overshooting is a consequence of our 
assumption that Q 2 can be treated as a large variable: in some terms 1/Q 2 serves as an expansion 
parameter. Such an approximation for these terms is invalid for small Q 2 : by analogy with the 
l/g 2 -terms, one may expect that, in the low-Q 2 region, one should substitute 1/Q 2 by something 
like 1/(Q 2 + m 2 p ). In fact, the change must be even more dramatic, since, in the Q 2 — > limit, only 
the anomaly term (i.e., the uncorrected contribution of the triangle diagram) should survive, while 
all other terms should vanish. Hence, using the limiting form 1/Q 2 , one appreciably overestimates 
the corrections for Q 2 ~ 1 GeV 2 , and this produces enlarged values for F 77 . 7r o(Q 2 ). 

It is also worth noting that the local duality approximation 



, so I 

*UF^AQ 2 ) = / P PT (s, Q\ q 2 = 0)ds= (8.4) 
'o l + Q z /s 



(applied to the original, uncorrected perturbative spectral density ( |3.6|) exactly reproduces the 
Brodsky-Lepage interpolation formula ( |2.24| ) if the duality interval Sq assumes the value Sq = 
47r 2 / 2 pa 0.67 GeV 2 , dictated by the local duality for the two-point function (see f50|j). 

In the region Q 2 > 3 GeV 2 , our curve for Q 2 F^* 7r o(Q 2 ) is practically constant, supporting 
the pQCD expectation ( [2.19|) . Comparing the absolute magnitude of our result with the pQCD 
formula, we conclude that it corresponds to the estimate / ~ 2.4 for the /-integral. Of course, 
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Figure 11: Combination Q 2 F jrY * 7r o(Q 2 )/4:iif n as calculated from the QCD sum rule (solid line), 
p-pole model (short-dashed line) and Brodsky-Lepage interpolation (long-dashed line). 

this value has some uncertainty: it will drift if we change our models for the bilocals (photon 
distribution amplitudes). The strongest sensitivity is to the choice of the photon distribution 
amplitude y? 7 (?/) in the non-diagonal correlator (|2.7|). However, it should be emphasized that even 
switching to the constant flat form 4> y (y) = 1 does not increase our result for I by more than 20%. 
The basic reason for such a stability is that the singular 1/q 2 factor from the relevant contribution 
in the original sum rule (|3.12 ) is substituted in ( j2.7|) by a rather small factor k/m 2 In fact, having 



precise data, one can obtain information on the shape of the photon distribution amplitude (p 7 (y) 



(ECU) 



Recalling that I as = 3 and I = 5, we conclude that our result / ~ 2.4 favours a pion 
distribution amplitude which is even narrower than the asymptotic form. In particular, if we 
parameterize the width of ip^(x) by a simple model if w (x) ~ [x(l — x)] n , we find that / = 2.4 
corresponds to n = 2.5. The second moment (£ 2 ) (£ is the relative fraction £ = x — x) for such a 
function is 0.125. This low value (recall that (£ 2 ) as = 0.2 while (£ 2 } cz = 0.43) agrees, however, 
with the lattice calculation [^TJ and also with the recent result obtained from the analysis of 
a non-diagonal correlator. 
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9 Conclusions 



Our basic goal in the present paper was to develop a regular QCD sum rule approach to the 
calculation of the transition form factor F 1J * 7r o(Q 2 ). Our results support the expectation that the 
(3 2 -dependence of the transition form factor F^* n o(Q 2 ) is rather close to a simple interpolation 
between its Q 2 = value (determined by the ABJ anomaly) and the large-Q 2 pQCD behaviour 
F{Q 2 ) ~ Q~ 2 . Moreover, the QCD sum rule approach enables us to calculate the absolute normal- 
ization of the Q~ 2 term. The value produced by the QCD sum rule is close to that corresponding 
to the asymptotic form (p^f(x) = 6/ 7r x(l— x) of the pion distribution amplitude. Though a detailed 
comparison with experimental data is beyond the scope of this paper, we would like to mention 
that our curve for F 77 . 7r o(Q 2 ) is in satisfactory agreement with the CELLO data || and in good 
agreement with the preliminary high-Q 2 results from CLEO II. We interpret our findings as a 
theoretical evidence that <p n (x) is a rather narrow function. 
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A Alpha-representation and asymptotic behaviour of the 
three-point function 

To study the asymptotic behaviour of the perturbative amplitudes in the limit when some mo- 
mentum invariants are large, one can use different types of integral representations for the relevant 
diagrams. For the purposes of a general analysis, one of the most effective approaches is that using 
the the "alpha-representation" for the relevant Feynman integral. To get the alpha-representation, 
one should write the denominator of each propagator of the Feynman diagram as 

1 f°° 

— ~ ~ = i exp{ia a (kl - m 2 + ie)}da a , (A.l) 

m z a — k* — le Jo 

where a numerates the lines of the diagram, and then take the resulting Gaussian integration over 
all the virtual momenta k a . As a result, for each diagram contributing to T(g 1; g 2 ) (see Fig.|l]a), 
one gets the expression having the following structure: 

\ \ ■ 2^o(«) . o^l( a ) • 2^2(a) ^ 2 1 

G{a,q 1 ,q 2 ;m (7 )exp<ip + iq 1 + iq 2 - i ^ a a {m a - le) \ , (A.2) 

{ D{a) D{a) D{a) V J 
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where d is the space-time dimension, P(c.c.) is the relevant product of the coupling constants, z is 
the number of loops of the diagram; D,Q,G are functions of the a-parameters uniquely determined 
by the structure of the diagram. In particular, D(a) is a sum of products of the a-parameters, 
with the number of the a-factors in each term of the sum being equal to z. In our case, all the 
functions Ai(a) are also the sums of products of the a-parameters, with z + 1 parameters a a in 
each product. Hence, D(a) and all Ai(a) are positive for positive a's. The preexponential factor 
G(a, qi, q 2 ; m a ) is a polynomial in a's, p 2 , q\ and q\. 

In the region where one of the momentum variables pj is large, all the contributions having 
a power-type behaviour on that variable can only come from the integration region where the 
relevant Ai/D factor vanishes: if A^/D is larger than some constant p in the region of integration, 
the resulting contribution is ~ exp(ip 2 p), i.e., it is exponentially suppressed. 

When all A^s are non-negative, there are two basic possibilities to arrange Ai/D = 0. In 
the first case, called the "short-distance regime", Ai vanishes faster than D when some of the 
a-parameters tend to zero (small a correspond to large virtualities k 2 , i.e., to short distances). 
The second possibility, called the "infrared regime", occurs if D goes to infinity faster than Aj 
when some of the a-parameters tend to infinity (large a correspond to small momenta k, i.e. 
to the infrared limit). One can also imagine a combined regime, when Ai/D = because some 
a-parameters vanish and some are infinite. 

In fact, there exists a simple rule using which one can easily find the lines a whose a-parameters 
must be set to zero and those whose a-parameters must be taken infinite to assure that Ai/D = 0. 
First, one should realize that Ai/D = means that the corresponding diagram of a scalar theory 
(in which G = 1) has no dependence on the momentum invariant p 2 . As the second step, one 
should incorporate the well-known analogy between the Feynman diagrams and electrical circuits 



53 1 : the a CT -parameters may be interpreted as the resistances of the corresponding lines a. In 



other words, a a = corresponds to the short-circuiting the line a while a CT = oo corresponds to 
its removal from the diagram. Hence, the problem is to find the sets of lines {a}sD, W}ir whose 
contraction into point (for {&}sd) or removal from the diagram (for {u}ir) produces the diagram 
which, in a scalar theory, does not depend on p 2 . 

Thus, the rule determining possible topological types of the short- distance factorizable contri- 
butions is the following: if the part of the diagram corresponding to a short- distance subprocess 
is contracted into point, the resulting effective diagram should have no dependence on the large 
momentum invariants. 

The simplest situation is when the short-distance part coincides with the whole diagram. This 
configuration is allowed for any relation between q\, q\ and p 2 . In this case, all the currents J^X), 
J u {0) and ja(Y) (see Fig.|[) are close to each other, i.e., all the intervals X 2 ,Y 2 , (X — Y) 2 are 
small. However, if, say, the variable q\ is small, the dependence on large variables q\ and p 2 can be 
eliminated by contracting into point a subgraph containing the vertices corresponding to momenta 
q 2 and p. In this situation, the interval Y 2 is small while X 2 and (X — Y) 2 are large. Such a 
configuration is sensitive to the long-distance effects in the ^-channel. They can be described 
by introducing distribution amplitudes for the gi-photon. Finally, another interesting situation 
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is when both the photon virtualities are much larger than p 2 . Then there exists a short-distance 
subprocess which includes only the photon vertices: the interval X 2 is small, while Y 2 and (X—Y) 2 
may be large. In this case, there are long-distance effects in the axial current channel, which are 
usually described/parameterized by the pion distribution amplitude. 

B Calculation of some useful momentum integrals 

We calculate integrals over the momentum using the dimensional regularization. The basic integral 
is well-known: 

I(L r) = f d D v {p2Y = '(-rV 4 - D r(r + D/2)r(Z-r-J/2) 
{ , ) J ap [p2 + S} L (4ir) D / 2 T(L) T(D/2) (S) L — D I 2 ' 1 ' 

where D = 4 — 2e and d D p = d D p/ (2n) D . In fact, it is more convenient to express the results in 
terms of the integral 

For our calculations, we need the following integrals: 

J F (p 2 ) a (p-q) 2 P Jo T{a)T{(3)J F [{p - q) 2 + S] L V ; 

where L = a + /3, q = qx, S = q 2 xx, x = 1 — x. Omitting, for a moment, the integration over x, 
let us concentrate on the integrals 

{J(L,n), J p (L,n), J pa (L,n), ...} = / d D p ^ Y ™^ P f^ L ■ ' } ■ (B.4) 

Shifting the integration variable, we expand them in a standard way: 

(p + qY) n = (qY) n + n{qY) n '\pY) + ^ ~ ^ (qY)^ 2 (pY) 2 + ■■■ . (B.5) 



Now, the integration over d D p is straightforward and, keeping the terms up to Y 2 , we obtain 

J(L, n) = (q.Y) n R(L, 0) + w(w ~ ^ (qY) n ~ 2 Y 2 R(L, 1) + 0(F 4 ), (B.6) 



+ nY p UqY)^R{L, 1) + {n ~ ^ ~ 2) {qY) n ^Y 2 R{L, 2) + 0(Y 4 ) ) , (B.7) 
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J pa (L,n) = q p qA(qY) n R(L,0) + 



^-^l(qY) n ~ 2 Y 2 R(L, 1) + 0(F 4 )) + 



(n - l)(n -2) 



+ n(g p Y CT + ^)j(^ n #(£,1 
+ g pa (qY) n R(L, 1) + w(w ~ 1} (qY) n ~ 2 R(L, 2) {2F p F ff + ^F 2 } + 



(qY) n - 6 Y 2 R(L } 2) + 0(Y 4 ) + 



+ 



n 



(n-l)(n-2)(n-3)_ B 



(gF)"- 4 J R(L, 3) {F P F CT F 2 + 0(F 4 )} . 



(B.8) 



C Equations of motion and p-meson distribution ampli- 
tudes 



Here we demonstrate how one can use equations of motion to obtain relations between the moments 



of the p-meson distribution amplitudes. A similar analysis was done in refs. J45| and |54|, |55 
Consider the identity: 

(0|^ a (0)(iV-m)^Vv(*> =o, (c.i; 
where iV p a = (id z + gA(z))p a . Applying the Fiertz transformation we rewrite ( p.l|) as 

VJz)- 



(z'V -m)p a S (z) + [(i V - m) 75 J P(z) + [(z V - m) 7^ 



/3a 



1 r 



(iV - m) 7^75 + - (iV - m) 07,5 



/3a 



T M(5 (z) = 0. 



(C.2) 



/3a ^ ' ' 2 

To obtain a relation between the distribution amplitudes, one should substitute in (CJ.2) the 
expressions for the bilocal matrix elements like (|b.8|),(]Sl^,( |6.3q ),( |b.3^ ),( |b.53| ), differentiate them 
with respect to z and take z 2 = 0. By contraction with [a vp ] a p y we pick out a combination of the 
V-, A- and T-projections. There are three independent tensor structures 

Zy^-p Z p E u , (^^) (ZvPp ZpPi/Jj Pu^-p Pp^-vi 

and, as a result, we get three systems of equations: 



f*m 2 / d(3(3n( [x 3 (3 + x 2 ] n ~ 1 ), A + 



+ V p A C AA n{x n - 1 ) AA - 2C VA f v p m p {x n ) VA + -f v p m p C [V2]i (*"> [y2]4 (C.3) 



f p A C A ^(x n ) M = a vl fJm p C [vl]4 (x n ) [vl]4 + -fpmpC [V2 ]4(x n+1 } [V2W 



(n + 2) (x n ) A = - f p m p (x n+1 ) v - a vi f p m p (x n ) vi 



f 3 A p £ dppn([x 3 p + x 2 ] n - 1 ) 3A 



+ fs P / d(3(3n{ [x 3 p + x 2 ] n - 1 ) 3V + Am q f (x n ) T 



(C.4) 



(C.5) 
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where we have used the notation: (x n ) A = Jq dxx 1l (p p i (x, /x 2 ), etc. It should be noted that 
eqs.( |C.3| ),( |C.4| ) are new. Eq. (|C.5| ) was derived also in pE5fl , but the constant ayi was missed 
there. In the chiral limit, we can neglect the last term in ( p.5| ). Taking the infinite limit for the 
renormalization parameter, fi 2 — > oo in ( JCJ.5[ ) , we obtain the equation relating the moments of the 
asymptotic leading-twist distribution amplitudes: 

(n + 2) ft (x n ) A = - f p m p (x n+1 ) v - a vl f p m p (x n ) vv (C.6) 



The asymptotic distribution amplitudes ipy (x), ^Pa( x ) were originally obtained in [f|8 



<p a v s = ^{l-2xx), p a A s = 6xx, (C.7) 

Taking into account the normalization conditions ( |6.10| ), we conclude that there exists only one 
solution: 

<p™ 1 = 60xx(2x-l), (C.8) 

with 

(/ A ) « = _^_£ j (ayir= (C.9) 



Note, that ^Pvi^vii gi ven by eqs. ( p.7|) and ( |C.8| ), obey the condition that, for a longitudinally 
polarized p°-meson (i.e., when e^ =0 ~ ip a /m p + O {m p /p z ) ,asp z — > oo), the leading-twist part 
in eq. (|6.8|) provides the well known asymptotic twist- 2 vector wave function (cf. eq. fl6.4j )). The 



value of f p was also calculated within the SR method [f4~8| , and the result is in a good agreement 
with that dictated by equations of motion. 
Substituting eq. ( |C.9|) into eq. (|C.6|) , we get 

^±^- (x n ) A = (x n+ \ + 1 (x n ) vl . (CIO) 



This formula is an analog of the well known relations [|56| (see also [ 55J, |57j| ) between the moments 
of the spin-dependent structure functions gi(x), g2{x). 

The asymptotic form for the next-to-leading two-body distribution amplitudes can be directly 
extracted from the corresponding correlators ( |5.2| ) : 

if A \ = 30y 2 y 2 , f[vi ]4 = 420(^) 2 (2y - 1), ^ 2]4 = 30y 2 y 2 . (C.ll) 
Using ( C.4 ) for n = and n = 1 gives the following relations: 

5 

C[V2]4 = — C*A4 and C[ V i]4 = ^ C*A4- (C.12) 

In the main text of the paper we gave the estimates ( |6.34j) for Ca<± and C[vi]a based on the local 
duality arguments. These estimates satisfy the second relation above. The first one can be used 
to estimate C[V2]4- 
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Combining the SR method for the correlators ( |5.2| ) and expansion of the relevant composite 
operators in the conformal basis[] one can calculate the corrections to ( |C.7|) - (|C.9| ). It should be 
noted that this expansion should be consistent with the equations of motion ( |C.3| ) - ( |C.5| ) . Such a 
program for the case of the nonleading 7r-meson distribution amplitudes was performed in HR K4 
and it was shown there that the deviation from the asymptotic form is small. 



D Some properties of traceless combinations 

To construct the orthogonal projection operators Pi n ) onto the subspace of traceless symmetric 
Lorentz tensors of rank n, we use the techniques similar to those in Here we list some useful 
formulas concerning these projectors as well as some contractions that appear in the paper. 



By definition, for an arbitrary Lorentz tensor T we have [f]|| : 



P(n)T 



It is straightforward to derive the formula 



p Hx—IMi r T V\...V n 
r iri)v x ...v n 1 



fll... )l„ 



in in 



(D.l) 



(D.2) 



where T ul ^ Mn ~ 1 a is now traceless and symmetric in its first n — 1 indices and [/^] means that the 
corresponding index is absent. Choosing rln-fn-ij-a = s a j^i _ _ _ g^ 11-1 } we have: 



n 



i=l 

1 n 



Making use of the Nachtmann's [58] contraction 

^...^"{gf ...qp} n = 



' 1 2 \ n / 2 

gi£ _ j 



Cn(^) 



(D.3) 



(D.4) 



and some recursion relations for the Gegenbauer polynomials C^(rj) fl59fl , one can derive the for- 
mula: 



1 9 



z" 1 . . . z""" 1 {q«q» 1 . . . gf 1 " 1 } = - — z" 1 . . . z*> {gf 1 . . . q^} 

1 

n 



n dz 

'z' 



(D.5) 



where ?7 = i (qiz) / ' \J —z 2 q\r = —iy—z 2 qf/2. 

Using ( p.3| ) - ( p.5[ ), one gets for an arbitrary 4- vector s: 



n 2 



(D.6) 



7 To one loop accuracy, only the operators with the same conformal spin mix under the renormalization. 
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E Contact terms 



Here we derive eg. ( |7.11 ). Before considering the relevant contraction, let us note that, incorpo- 
rating the relation 



{zd} n = (zd){zd} n ' 1 - { 1 {zd} n ~\ 

4 n 

and neglecting higher twist contributions, one can substitute the original correlator (fT. 



(e.i; 



by 



As a result, 



n At { m ...^ n }(gi) = n MAtl { At2 ... At „}(gi) H . 

n M { Atl ... Mn }(?i) 9 W i - n W { /12 ... M „}(?i) 

' eke"^ (0|T{J M (x) u(0)M^ . . . a Mn }n(0)}|0) 



-- / rfa;e-^(0|T{J M (a:)n(0)57MR 2 ---^}^(0)}|0) 



+- / dxe~^ x (0|T{J M (x) M (0) 7 ,R 2 ...^J^(0)}|0) 



+-g le / dxe-^(0|T{J M (x)w(0) 7M 7 £ R 2 ...^ n K0)}|0), 



where we have made use of the identity: 



dxe- iqiQ 



(0|T{J m (x)m(0)«9 £ iM0)}|0) + (0|T{J M (x)M(0)r,9 £ «(0)}|0> 
= iq u J dxe~ iqiX (0|T{J M (x)M(0)r\i(0)}|0). 
Applying ( [7.3| ) to ( |E.3| ) and ( |E.4| ), and integrating by parts in ( |E.3| ) we get 

(E.3) = 2Hr- 1 {gi, 2 ...gi M , l }(^). 

Now, taking into account that (0\u(0){d^ 2 . . ■ d^ n }u(0)}\0) = for all n, we obtain: 

(E.4) = . 



(E.2) 



(E.3) 
(E.4) 
(E.5) 



(E.6) 



(E.7) 



(E.8) 
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